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ABSTRACT 


A  practically  important  regularization  of  the  Navier-Stokes  equations  have  been  ana¬ 
lyzed.  As  a  continuation  of  the  previous  work,  we  study  in  this  paper  the  structure  of  the 
attractors  characterizing  the  solutions.  Local  as  well  as  global  invariant  manifolds  have 


been  found.  Regularity  properties  of  these  manifolds  are  analyzed. 
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Chapter  1 
Introduction 


In  [15]  we  have  presented  an  analysis  of  a  practically  important  form  of  the  regularized 
Navier-Stokes  equations.  We  also  presented  a  theorem  for  the  convergence  of  the  solution  of 
such  system  to  the  solution  of  the  conventional  Navier-Stokes  system  as  the  regularization 
parameter  approaches  zero. 

In  this  paper  we  will  present  an  analysis  of  the  structure  of  the  attractors  associated 
with  the  regularized  system.  We  will  present  in  particular  theorems  for  stable  and  un¬ 
stable  manifolds  associated  with  each  periodic  solution  and  establish  their  analyticity  and 
invariance  properties.  The  main  machinery  needed  for  these  invariant  manifold  theorems 
are  the  analyticity  properties  of  the  nonlinear  semigroup  and  its  Frechet  derivative  and 
spectral  theorems  for  the  monodromy  operator.  These  results  are  established  in  section  3. 

In  the  section  4  of  the  paper  we  establish  the  existence  of  a  global  attractor  for  the 
system  and  prove  its  compactness.  We  also  note  certain  bounds  on  the  attractor  which  are 
uniform  on  the  size  of  the  regularization  parameter.  We  then  prove  the  existence  of  global 
(inertial)  invariant  varieties  containing  this  attractor.  Such  a  global  invariant  manifold 
theory  is  proposed  in  [6]  for  certain  class  of  semilinear  evolution  equations.  Motivation 
for  such  study  of  course  comes  from  the  famous  paper  of  E.  Hopf  [8].  We  then  study  the 
regularity  of  the  inertial  manifolds  and  obtain  sufficient  condition  for  them  to  be  Cl. 
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Chapter  2 

Governing  Equations  and  Functional 
Framework 

In  this  chapter  we  will  briefly  outline  the  mathematical  framework  used  in  this  paper. 
For  detail  proofs  of  the  relevant  theorems  see  [15].  We  regularize  the  conventional  Navier- 
Stokes  equations  by  adding  a  fourth  order  operator  (Laplacian  square)  with  an  artificial 
dissipation  parameter  e.  In  addition  to  the  prescribed  initial  field  u0  and  Dirichlet  bound¬ 
ary  condition,  we  also  prescribe  the  Laplacian  of  the  velocity  field  at  the  boundary  to  be 
zero.  Let  fl  C  Rn,n  <  6  be  a  bounded  open  set  of  class  Cr,r  >  4.  The  problem  is  to  find 
(«,  p)  :  fix  (0,  oo)  —yRnxR  such  that 


du 

— — |-  cA2u  —  i/Au  +  (u  •  V)u  +  Vp  =  /  , 
at 

in  f2  x  (0,  oo)  , 

(2.1) 

< 

a 

II 

o 

in  f2  x  (0,  oo)  , 

(2.2) 

u|an  =  0,  Au|an  =  0, 

(2.3) 

«(x,0)  =  tt0  , 

in  fl . 

(2.4) 

Here  u  >  0  is  the  coefficient  of  the  kinematic  viscosity  of  the  fluid  and  /  is  a  prescribed 
vectorfield. 

Let  us  introduce  the  following  function  spaces: 

j( n)  =  {u  :  fl  ->  Rn\  u  e  C*( n),  u|an  =  0,  Au|an  =  0,  div  u  =  0}, 

H  =  {«  :  fl  — >  Rn ;  u  G  L2( fl),  div  u  =  0  u  •  n|an  =  0}, 

V  =  {u  :  n  ->  Rn,  v  e  H2{ n),  div  u  =  0,  u|an  =  0}. 
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[ 


Here  we  denote  by  Hm{Q),  the  Hilbertian  Sobolev  space  of  (square  integrable)  vector- 
fields  whose  distributional  derivatives  up  to  order  m  are  square  integrable.  These  spaces 
are  endowed  with  the  inner  product 

(«,v)tfm(n)  =  (-Da“,-Dav)L2(  n) 


and  the  norm 


llulltfm(n)  -  (  \'Dau\\2Li(n j)1?2  . 

|a|<m 

It  ran  be  shown  [15]  that  the  spaces  H  and  V  are  respectively  the  completion  of  j( fl)  in 
the  norm  L2( fl)  and  H2(fl). 

The  space  H  is  endowed  with  the  inner  product  (u,  t>)x,2  and  norm  |«]  =  ( u,« )z'a  .  One 
can  easily  verify  that  the  norm  induced  by  H2[fl)  and  the  norm  ||Au||£,2(n)  are  equivalent 
in  V.  We  then  denote  ||u||  =  ||A«|(x,2(n)  =  (u,u)]/2  as  the  norm  in  V  derived  from  the 
inner  product 

,  ^  _  vv  d*u  d*v  ^ 

f^dxxdxi'  dxidx) ' 

Let  us  now  characterize  a  linear  self  adjoint  positive  operator  A  (which  we  call  the 
dissipation  operator)  using  the  following  fundamental  linear  problem  :  Find  (u,  p)  :  fl  — ► 
Rn  x  R  such  that 


A2u  +  Vp  =  /  , 
Vu  =  0, 


in  fl , 
in  fl , 


(2.5) 


[  «|an  =  0  ,  Au]an  —  0  . 

This  equation  characterizes  a  linear  operator  relating  u  to  /  and  is  defined  on  j(O).  Let 
us  denote  the  Friedrich’s  extension  of  this  operator  as  A.  This  operator  can  be  defined  as 
follows:  We  define  a  positive  definite,  V-elliptic  symmetric  bilinear  form  a(-,  •)  :  V  x  V  — *  R 

by 


«(«>»)  =  E( 


d2u  d2v 


-)■ 


dxidxi'  dxxdxi 

Then  by  Lax-Milgram  lemma  we  obtain  an  isometry  A  G  C(V\V')  as 
<  Au,  v  >v'xv~  a(u,v)  =<  f,v  >v'xv,  V  v  6  V, 
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and  Au  —  /  G  V'  =  t{V\R).  We  then  define  D(A)  as  follows:  for  f  (E  H  c  V',  3  u  G  V 
such  that 

a(u,»)  =  V»e7. 

We  then  denote  «  G  £>(^4).  We  thus  have  A  G  Z(D{A)\H)  fl  £(V;V').  The  operator  A 
defined  above  is  closed  with  D[A)  dense  in  V  <Z  H.  From  this  it  is  easy  to  conclude  that 
A  is  self  adjoint.  Since  the  continuous  form  a(-,  ■)  is  positive  definite  we  deduce  from  a 
theorem  of  Lions  [13]  that  D(A l!2)  =  V  and 

a(u,v)  =  (A1/2ti,A1/2v)  V  u,  v  e  V  . 

This  implies, 

||  rx|| 2  =  a(u,  u)  =  \Al!2u\2  V«G  F. 

We  have  in  fact  A  =  A\A2  :  u  —*  f  so  that  A  is  an  isomorphism  from  D(A)  onto  H 
with  Al  the  Stokes  operator  and  A2  the  Friedrich’s  extension  of  the  Laplacian  operator[l5]. 
This  gives 

u  G  D{A)  =  {«  G  Lf4(f2);«|an  =  0,  Au|an  =  0,  divu  =  0}. 

We  also  have  the  following  estimate  for  solution  (u,  p) 

||«||/r«(n)  +  ||p||ffi(n)/K  <  coll/ll h  ■ 

In  consequence  of  the  relationship  Au  =  /,  3/?i,/?2  G  R+  such  that 

A||tt||jr«(n)  <  \Au\  <  &Nlir«(n)  ,  VuG  £>(,4).  (2.6) 

By  Rellich’s  Lemma  [1]  A~l  as  a  mapping  in  V'  (or  H)  is  compact.  Hence  the  spectrum 
of  operator  A  consists  of  real  eigenvalues  p;-  of  finite  multiplicities  and  can  be  ordered  as 

0  <  Pi  <  P2  <  •  •  ’  ;  Pj  — ►  +oo  as  j  —*  +oo, 

with  accumulation  possible  only  at  infinity.  The  self  adjoint  operator  A  possesses  an 
orthonormal  set  of  eigenfunctions  complete  in  V  (or  H). 

A<f>j  =  p,0y  ,  <(>,  G  V (or  Z?(i4)),  V  j  .  (2.7) 
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If  we  identify  H  with  its  dual  H'  using  the  Riesz  representation  theorem  then  we  get 
the  following  continuous  dense  and  compact  embedding  structure  : 

D[A)  CVCH  =  H'CV'  C  D(A)'. 

Let  us  now  define  the  trilinear  form  associated  with  the  inertia  terms: 

n  -  q 

b(u,v,w)  =  £  /  ViDiVjWjdx,  D,  =  — — . 

«,j= iJn  dx* 

It  can  be  easily  shown  by  applying  the  divergence  theorem  and  noting  that  v  and  to  have 
zero  trace, 

b(u,v,w)  =  — 6(u,to,v)  ,  Vu,t),t«GV 

and 

6(u,  v,  v)  =  0  ,  V  u,  v  G  V  . 

By  the  application  of  Holder  inequality  and  Sobolev  embedding  theorem  we  can  show 
that  6( •,•,•)  is  trilinear  continuous  on  Hmi( O)  x  i/m2+1(fi)  x  iLm3(n),  m,  >  0: 

|6(«,v,  to) |  <  c0||u||fl'mi(n)|!»||^2+i(n)l!«J||w-«3(n)  , 


mi  A  m2  +  m3  >  -  if 
n  .. 

mx  +  m2  +  m3  >  —  if 


i  =  1,  2, 3  and 
for  some  i . 


(2.8) 


In  particular,  for  n  <  6,  6  is  a  trilinear  continuous  form  onV’xV’xV.  When  O  is  bounded, 
the  following  interpolation  inequality  holds  [14]: 


|  ti  jj(1  —  9)mj  +0ni2  ^  C 


Vu£  Hm2( fl),  mi  <  m2,  6  G]0, 1 


From  the  above  we  will  derive  in  particular, 

|6(u,v,to)|  <  cx|«|1/2||tt|j1,,!iit;!j1',2jJ4t;|1/2|tt7|,  Vu  G  Vr,o  G  £>(A),to  G  H  (2.10) 

|6(«,  w, «>)[  <  c2|«l1/2IMI1/2M1/2IMI1/2IHI>  v  u,v,w  g  v  ,  (2.11) 
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where  C!,c2  are  positive  constants.  By  virtue  of  (2.8),  we  know  that  the  above  inequalities 
are  valid  for  space  dimension  n  <  6.  The  estimate  (2.8)  enables  us  to  define  (using  Riesz 
representation  theorem)  a  bilinear  continuous  operator  B  from  Hm'(Q )  x  Hm2  +  l(fl)  into 
(//m3(f]))'.  In  particular,  for  u,v,w  E  V,  B(u,v )  E  V'  will  be  defined  by 

<  B(u,v),w  >v<xv~  b(u,v,tu)  ,  V  w  E  V  .  (2.12) 

Let  us  note  that  a  linear  operator  Ai  from  V  onto  H  can  be  defined  as 

(AxU,v)  =  (Vu,Vv),  V  u,  v  E  V  (2.13) 

and  A\U  =  —PHAu,  Vu  E  V.  In  fact  Ax  is  the  Stokes  operator  associated  with  the 
conventional  Navier-Stokes  equations.  PH  is  the  orthogonal  projector  in  L2(fl)  onto  //. 
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Chapter  3 

Local  Invariant  Manifolds 


In  this  chapter  we  will  establish  the  hyperbolicity  of  periodic  solutions.  Existence 
theorem  for  the  periodic  solutions  is  provided  in  [15].  Here  we  will  study  the  orbits  nearby 
each  periodic  orbit.  We  will  prove  in  particular  the  existence,  uniqueness  and  analyticity  of 
stable  and  unstable  manifolds.  Such  results  for  conventional  Navier-Stokes  equations  have 
been  proven  in  [20].  Let  us  consider  a  perturbation  about  a  general  time  dependent  smooth 
and  bounded  field  ( U[x ,  t ),  P{x,  t)).  We  assume  that  this  basic  field  satisfies  the  governing 
equations  and  the  boundary  conditions.  Let  us  introduce  u  =  U  +  v  and  p  =  P  +  q  in 

(2.l)-(2.4).  Then  the  perturbation  solution  ( v,q )  satisfies. 

(  dt) 

I  —  4-  cA2»  -  uAv  +  ( U  •  V)tr  -j-  (v  •  V)U  +  (v  ■  V)v  +  Vq  =  0, 


V  •  v  =  0,  in  fl  x  (0,  oo), 

«|an  =  0  Aw|an  =  0, 


(3.1) 


(  v(x,0)  =  Vq,  in  fl. 

Let  us  now  rewrite  this  system  as  an  equation  of  evolution  in  the  Hilbert  space  V .  This 
can  be  achieved  by  applying  the  projection  operator  Ph  on  the  system  (3.1).  Noting  that 
PH(Vq)  =  0,  we  get  the  evolution  equation  for  v: 

(  dv 

—  -r  eAv  -t  uA\V  +  Lu(t)v  +  B(v,  u)  =  0,  t  >  0, 

{  dt  (3.2) 


(  t?(o)  =  v0  e  v. 


Here  A  6  £(V;V')  is  the  dissipation  operator  defined  earlier  and  Ai  £  £(V;IJ)  is  the 
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Stokes  operator.  We  can  use  the  Riesz  representation  theorem  to  characterize  Lu[t)  and 
B(-,  ■)  as 

<  Lv(t)v.  w  >ylyV=  b(U(t),v,w)  +  b(v,  U(t),  w),  V  w  €  V. 


and 


<  B(v,  v),  w  >v'*v=  b(v,  v,  t/j),  X/weV. 


3.1  The  Cauchy  Problem  and  Associated  Semigroup 

In  this  section  we  will  derive  few  useful  properties  of  the  semigroup  generated  by  the 
dissipation  operator  -eA.  These  estimates  will  be  used  to  resolve  the  nonlinear  semigroup 
associated  with  the  regularized  Navier-Stokes  equations  as  well  as  to  establish  certain 
results  concerning  the  invariant  manifolds. 

Let  us  consider  the  Cauchy  problem: 

Problem  1  Find  v  6  C([0,  oo);  V)  n  (^(O,  oo;  V')  such  that 

dv 

— —  +  eAv  =  0  ,  t  >  0  , 

<  dt 

.  v(0)  =  v0  €  V. 

Before  studying  the  semigroup  S(t)  :  v0  —*  v(t)  associated  to  the  above  problem,  we  will 
document  certain  relevant  properties  of  the  resolvent  of  —  eA. 

Lemma  3.1  Let  B\  and  B2  be  two  Hilbert  spaces  defined  below  such  that  the  embedding 
B i  C  B2  is  continuous  and  dense.  Then  the  resolvent  of  the  dissipation  operator  -eA 
satisfies  : 

(i)  j|7?(A; -£>1)||£(B2.j52)  <  —  ,  for  A/  /  0,  where  Xr  =  ImX,  (3.3) 

(ii)  \\R(X-, <  T - p - ,  for  XR>  -epi  ,  (3.4) 

*R  +  £Mi 
M 

(iii)  \\R(X] -cA)\\c(b2,b,)  <  rr— - r ,  for  X  <E  E6,  X  -ea  ,  (3.5) 

A  -t-  (Cl 

where  £$  =  {A;  \arg(X  +  ea)  |  <  n/2  +  6,0  <  6  <  ir  j  2}  and  0  <  a  <  px ,  px  is  the  smallest 
eigenvalue  of  the  operator  A.  Here  we  will  take  the  spaces  B\  and  B2  as  either 
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(a)  By  =  V,  fl2  =  V" 


(b)  By  =  D(A),  B2  =  H. 

Proof:  (i)  let  us  write  u  in  By  as  an  expansion  of  the  eigenfunctions  of  A: 


U  =  22(u,<f>k)B2<f>k  ■ 


Since  <f>k  are  orthonormal, 


,  D/,  1 1 2  _  V'  I  :)  B7  1  -  j;  1 1  1 1 2 

,R(A;-eA)u !!«,-£  |A  +  e^|2  -  jA/|2  :i  llB2  ' 


Here  we  used  the  fact  that 


- - <  -  V  k,  A  =  Ar  4*  iA/ 

(Ar  +  ^)2  +  |A/I2  -  |A/|2’ 


Thus,  the  resolvent  of  -eA  in  B2  satisfies 


||i?(A; -cA)||£(b3;b2)  <  j^Ti  ,  for  A/ ^  0 


(ii)  Similarly,  it  follows  from 


_ \ _  <  1  V& 

(Ar  +  €M*)2  +  | A/|2  ~  (Ar  +  e^i)2 


that  for  Ar  >  —  efj.y, 


||^(A; -Ci4)||£(B2;B2)  <  t— 

aR 


+  CMi 


(iii)  Combining  (i)  and  (ii)  we  get  for  0  <  a  <  Hy 


|i/?(A;  — eA)||£(B2iB5)  <  |A  +  ea|  >  Wlth  ^R>  ea" 


From  this  we  show  by  analytic  continuation  that 


1  y/ 

||f?(A;  -iA)\\c(b2-b3)  -  |A  +  fa“  ’  for  A  G  X^~ea- 
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Theorem  3.1  —lA  generates  a  C°  semigroup  S(t)  such  that: 


(i)  l|S«llr(B,;Bl|  <  e-*“  ,  (>0, 

O'1)  l|5(0IU(Ba;fll)  <  je-£af,  t>  0, 

where  0  <  a  <  pi  and  c  depends  on  e.  Moreover  S(t )  can  he  extended  as  a  holomorphic 

semigroup  in  the  sector  A s  —  {z  ■  |ar<7.z|  <  6,6  =  tan-1  —,Rez  >  0}. 

ec 

Proof:  (i)  The  spectrum  of  the  operator  —eA  lies  on  the  negative  real  axis,  thus  the 
resolvent  set  p(— eA)  will  contain  the  positive  real  axis  and  from  the  estimate  on  the 
resolvent  in  Lemma  3.1,  we  have 

||JR(A; -£A)\\c{B2.B2)  <  - — -  for  XR  >  -ea,  a  <  ^  . 

A  r-\-  ea 

Hence  by  the  Hille-Yosida  theorem  [16],  —  eA  generates  a  strongly  continuous  semigroup 
S(')  in  Bi  and 

||5(t)|U(B2;B2)<e”to<,  t>  0 

(ii)  Using  the  estimate  (iii)  of  Lemma  3.1,  We  can  represent  S(t)  as  an  integral 

S{t)  =  -L  [  eXtR{\]  -eA)  dX 
2m  Jr 


where  T  is  a  smooth  curve  in  Y,s  consisting  of  two  rays  xe'e  and  xe~'e  ,  0  <  x  <  oo  and 
7t/2  <  6  <  7T.  T  is  oriented  so  that  A /  increasing  along  I\  Differentiating  the  integral  with 
respect  to  t,  we  get 

S\t)  =  —  f  XextR(X\ -eA)  dX  . 

2nt  Jr 

From  (3.5),  for  A  ^  —  ea  and  t  >  0 


|‘S,'(t)||£(B2;B2)  < 


Consequently, 


ell5(0ll£(B2;B,)  -  e||i45(t)[|£(B2;B2)  =  ||,5,(t)||£(B2;B2)  <  -«  eat ,  for  t  >  0. 
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Now  due  to  the  differentiability  of  S(t)  for  t  >  0  and  =  (  — eyl5(— ))n  =  (S'(  —  ))n, 

n  n 

we  have 

To  extend  the  semigroup  S' (2)  in  some  sector,  we  consider  power  series  for  S(z)  around  t 

=  s(t)  +  £  - 1)“ . 

n=  1  U‘ 

This  series  will  converge  in  Z{Bi\ B2)  for 

\z  —  t\  <  k  ^ j  ,  for  k  <  1,  t  >  0  . 

Hence  S(t)  can  be  extended  to  a  holomorphic  semigroup  S(z)  in  the  sector  A$  =  {2  : 

\arg  z\  <  6, 6  =  tan-1  — ,Rez  >  0}. 

ce 

□ 


Using  the  estimate  (iii)  of  Lemma  3.1,  we  can  define  the  positive  as  well  as  negative 
powers  of  A.  Since  A~x  is  compact,  the  spectral  resolution  of  the  self  adjoint  operator  A 
can  be  used  to  define  its  fractional  powers  in  a  simple  way.  We  thus  write,  Va  >  0 

OO 

A *«  =  £>««,  *,)*„  Vue  D(Aa), 

k=  1 

with  D(Aa )  =  {u  €  H  such  that  X^i  ^“Kti,  4>k) |2  <  00}. 


Lemma  3.2  For  a  >  0  and  t  >  0,  the  bounded  operator  AaS(t )  satisfy  the  estimates: 

II  £(B2;B3)  < 


M  ,  a 

—  for  0  <  t  <  - , 

t°  J  ep. 


a 


for  t  >  - , 

epi 


where  M  =  ( — )“. 

ee 

Proof:  Note  that  any  element  u  in  Bx  may  be  represented  as 


ti 


OO 

Y^(U^k)B74>k  ■ 


k=  1 
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Hence  we  may  write 


A*S(<)«  =  £ 

k~l 

Thus, 

||^a‘?(0IU(B2;B2)  <  max^e'^')  . 

CL 

Note  that  the  right  hand  side  of  inequality  admits  maximum  when  —  — •  This 

et 

proves  the  lemma  since  ^  >  0  is  the  smallest  eigenvalue. 


□ 


Let  PN  be  an  orthogonal  projector  in  V'  onto  the  finite  dimensional  subspace  of 
span{^x,  ■  •  • ,  <f>N}  and  Qn  =  I  —  Pn •  Note  that  Pn  and  Q n  commute  with  Aa.  The 
following  lemma  is  useful  in  a  later  section  and  can  be  proved  in  a  similar  way  to  the 
previous  lemma. 


Lemma  3.3  For  a  >  0  and  t  <  0,  the  bounded  operator  AaQNS{—t )  satisfy  the  estimates: 


|y40,Q^5'(-t)j|£(B2;B2)  < 


M 

H)“ 


a 


for - <  t  <  0  , 


£^n+i 


flN+lef>lN+lt  for  ~°°  <  t  < 


CL 


t^N+1 


OL  , 


where  M  =  ( — )“. 

ee 


□ 


Since  semigroup  S(t )  is  holomorphic  we  have  for  t2  >  >  0 

S{t2)  -  S{U)  =  [h  d~^-dt  =  f\-eA)S{t)dt. 

Jt  i  at  Jti 

This  gives 

||5(t2)  -  ‘S'(*i)||£(/T;V)  <  c  [  \\AS(t)\\c(H.V)dt  =  e  f  \\A3/2S(t)\\C(H.H)dt. 

Jti  J  t\ 

By  applying  Lemma  3.2  with  a  =  3/2  and  B2  =  H,  we  obtain 

!!*s(t2)  -  5(«,)I|£(H;V)  <  Cli-T/i  -  -Tji)  ,  h  >  ty  >  0.  (3.8) 
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3.2  The  Characterization  of  the  Monodromy  Opera¬ 
tor 

We  will  now  characterize  the  evolution  operator  Z (■,•)  associated  with  the  Cauchy’s 
problem  obtained  by  linearizing  the  regularized  Navier-Stokes  equations  about  a  smooth 
time  dependent  basic  field.  Let  us  consider  the 

Problem  2  Find  v  G  C ([0,  oo);  V)  fi  C1((0,  oo);  V')  such  that 

dv  ,  .  ,  . 

- f-  eAv  +  uAiV  +  Lu(t)v  —  0,  t  >  0  ,  (3.9) 

dt 

v(0)  =  v0  e  V  . 

We  will  show  that  the  Problem  2  is  equivalent  to  the  following  integral  representation 
for  v(t): 

Problem  3  Find  v  6  C([0,  oo);P)  such  that, 

v(t)  =  S(t)v0  -  [  S{t  -  t)[i/Ai  +  LV[T)]v(T)dT  ,  (3.10) 

Jo 

v(0)  =  »o  G  V  . 

Theorem  3.2  Let  the  basic  field  satisfies  U  €  C([0,  oo);  H1^))  and  be  bounded.  Then 
the  Problem  S  resolves  the  Problem  2. 

Note  that  the  Stokes  operator  A\  is  a  linear  continuous  operator  from  V  onto  H .  Further¬ 
more,  the  linear  operator  Lu(t)  is  characterized  by 

<  Lv{t)v,w  >v'xV=  b{U(t),v,xv)  +  b(v,  U(t),  to),  Vw  G  V. 

By  virtue  of  (2.8)  we  obtain  the  following  lemma. 

Lemma  3.4  IfUeH^Q),  LveL{V-,H). 
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Proof:  Recalling  the  estimate  for  6(-,-,-)  in  (2.8),  we  take  mi  =  2  and  m2  =  m3  —  0  to 
get 

|6(w,  U,w)\  <  Cillvll^ijnjlir/Uff^fijllt^llijjn)  ,  V  t u  £  H  . 

A  similar  estimate  holds  for  b(U,v,w).  Thus 

\(Lirv,w)Li(n)\  <  (ci  +  c2)||t;||v||^||//*(n)||«JIU2(n)  ,  V w  G  H  . 

Here  we  note  that  the  norm  induced  by  H2(Cl)  is  equivalent  to  the  norm  in  V.  Setting 
w  =  Lyv  <E  H  we  get, 

\\Luv\\li{(i)  <  (Cl  +  c2)||C/||Wi(n)||v||v- 

□ 

In  order  to  prove  Theorem  3.2,  we  need  first  to  establish  certain  properties  of  linear 
operator  K  defined  as 

[Kr](0  =  -  [  S(t  -  r){i/Ai  +  Lv{t)]v{t)(It  . 
jo 

Lemma  3.5  Let  K  be  the  linear  operator  defined  above.  Then  for  sufficiently  small  Tx, 
K  is  a  contraction  in  C([0,  Ti);  V) .  Moreover,  K  can  be  extended  as  a  contraction  in  the 
Banach  space  B  defined  as:  B  =  {t  — »  »(*);»(£)  continuous  in  V  for  t  e]0, Ti]  and 
bounded  in  V  },  with  norm 

IMIb  =  sup  ||t1/2u(t)||v  . 

<e(o,r,) 

□ 

Lemma  3.6  For  v  £  C([0,  Ti);  V),  the  time  derivative  [£»]'(•)  £  C((0,  Tj);  V').  Moreover, 
when  v  €  B  the  map  t  — *  [X’w]'(t)  is  continuous  from  ]0,  Ti)  into  V'. 

The  proofs  for  the  above  two  lemmas  are  similar  to  those  for  the  conventional  Navier- 
Stokes  equations  [20].  Note  that  [X’vJ(t)  €  C([0,Ti);7)  and  it  has  continuous  right 
derivative  [Kt>]'+(i)  €  C((0,  Ti);  V').  This  implies  [Kw](f)  is  strongly  differentiable  and 
[X’wj'(t)  =  [Kv]'+(i)  for  t  E  (0,  Ti)  (see  Zaidman  [23]). 
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□ 

Proof  of  Theorem  3.2:  We  have  shown  that  the  strong  derivative  [JCt>]'(<)  exists  and 
[X =  —eA[Kv\(t)  —  \uAi  +  Lu(t)}v{t).  This  gives 

(X»]'(t)  +  eA[Xv](£)  +  [uAi  +  Lu(t)}v(t)  =  0, 

with  [ XT t;] ( 0)  =  0  which  implies  the  representation  (3.10)  in  Problem  3  satisfies  the  differ¬ 
ential  equation  (3.9)  in  Problem  2.  Moreover,  due  to  the  properties  of  the  linear  semigroup 
S  ( t )  established  in  previous  section  we  have, 

5(t)»0  E  C([0,oo);  V)  n  ^((Ojoc);  V') . 

From  Lemma  3.5  and  3.6,  we  can  conclude  that 

T  S(t  -  t)[uAi  +  Lu(r)}v(T)dT  e  C({ 0,  T,);  V)  n  C\{0,  T V'). 

Jo 

Let  us  now  characterize  the  evolution  operator  associated  to  the  linear  differential 
equation  in  Problem  2.  From  the  definition  of  linear  operator  K ,  we  can  rewrite  (3.10)  in 
Problem  3  as 

[(/-  K)t;)j(t)  =  S(t)v0. 

Then  for  small  enough  Tu  the  operator  (I  —  K)  is  invertable  in  C([0,  Tj);  V)  and  in  B  since 
X"  is  a  contraction  in  these  spaces.  We  obtain  a  convergent  series  in  C([0, Ti);l^)  as 

»(■)  =  ElK-SlOo-  (3.11) 

n=0 

Hence,  the  solution  of  Problem  3  can  be  denoted  by  v(t)  =  Z(t,0)vo.  We  call  Z(t,  0)  the 
evolution  operator.  Note  that  the  convergence  of  this  series  ensures  the  uniqueness  of  the 
solution  to  Problem  3  (and  hence  to  the  Problem  2.)  Let  us  now  study  the  evolution  for 
t  >  t  by  prescribing  the  initial  data  at  t  =  r  in  Problems  2  and  3.  The  evolution  operator 
obtained(as  a  series)  in  this  manner  is  denoted  Z(t,r)  with  t  —  t  <  T\.  Here  T\  is  taken 
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small  enough  to  ensure  the  convergence  of  the  series.  Let  us  consider  for  0  <  r?  <  t  <  t, 
v(t )  =  Z(t,T)v(r)  with  »(r)  =  Z(t,t])v0.  That  is 

V{t)  =  Z{t,T)Z(T,T))v0. 


Due  to  the  uniqueness  of  the  solution  we  have 


»(*)  =  Z(t,r})v  o  =  Z(t,r)Z(r,  rj)v0. 


That  is 


Z(t,r j)  =  Z(t,r)Z(T,r ?),  0  <  r]  <  r  <  t . 

Iterating  this  kind  of  arguements  we  can  extend  the  definition  of  Z(t^,  tx)  to  t2  —  ti  £  [0,  oo). 


□ 


The  next  series  of  lemmas  provide  useful  regularity  and  compactness  properties  of 
evolution  operator  Z These  lemmas  can  be  proved  using  the  same  methods  used  in 
the  context  of  conventional  Navier-Stokes  equations  [20].  Note  that  here  the  initial  data 
is  prescribed  in  V  for  the  evolution  problem.  This  means  we  need  to  characterize  Z(-,-) 
as  an  element  in  £(V-,V).  Moreover,  the  bilinear  operator  characterizing  the  inertia  term 
has  the  property  B(-,  )  £  C{V  x  V;H).  Hence  we  need  to  extend  Z(-,-)  as  an  element  in 
C{H-V). 

Lemma  3.7  For  0  <  t  <  t  the  evolution  operator  Z(t,r)  satisfies  the  following  estimates: 

(i)  \\Z(t,T)\\c(v.v)  <  c^e01^  r)  , 

00  <  C,{1  +  ■ 

Here  c3,c4  >  0  and  <7i,<72  >  0. 


□ 
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Lemma  3.8  For  0  <  t2  <  t2  <  oo  we  have  as  t2  — ►  tf  ,Z{ti,tf)  — ►  I  strongly  in  £(V;V). 

For  t  T  the  TYlCLp  t  ^  ^  j  T  j  15  CO TXt  1  YXXt OXJLS  1 YX  the  \XYX\f Of*  171  O pc Y OLtOY  to po/o^y  O ^  ^  |  j  j  f~1 

i2(i?;  V^)  and  /oy  t  <  t  the  map  r  — ►  Z{t,r)  is  continuous  in  the  uniform  operator  topology 
of  £{V;V)n  £{H;V). 

□ 

Lemma  3.9  For  0  <  r  <  t  the  operator  Z(t,r)  6  C{V ;  V)  n  Z{H\  V)  is  compact. 

□ 

Note  that  the  evolution  operator  Z(t2,ti)  :  H  — >  V  is  compact  and  hence  the  spectrum 
of  this  operator  is  discrete  with  finite  multiplicity  and  accumulation  possible  only  at  the 
origin.  Moreover,  they  are  the  same  in  H  and  V . 

Let  us  now  specialize  our  study  to  the  case  where  the  basic  field  U  is  T-periodic  in 
time. 

Lemma  3.10  Let  the  basic  field  be  T-periodic  in  time.  Then 

(i)  Z(nT,0)  =  Z(T,0)n  ,Vn>  1, 

(ii)  The  spectrum  of  Z(T  +  t0,t0)  is  independent  of  t0  >  0. 

□ 

We  will  call  the  operator  Z(T,  0)  the  Monodromy  operator. 

3.3  The  Nonlinear  Semigroup 

In  this  section  we  will  characterize  nonlinear  semigroup  associated  with  the  nonlinear 
syatem  (3.1).  We  will  define  in  particular  the  time  T-map  which  relate  the  initial  data  v0 
to  the  solution  v  at  time  T.  In  addition,  we  will  prove  that  the  uependence  of  v  in  the 
initial  data  is  Frechet  analytic.  This  last  result  will  be  used  in  next  section  to  establish 
the  analyticity  of  the  local  invariant  manifolds.  Let  us  now  consider  the  evolution  form  of 
the  regularized  Navier-Stokes  system. 
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Problem  4  Find  v  €  C^O,  oo);  V)  n  C1(0,c»;  V')  such  that 

dv 

—  +  eAv  +  uAyV  +  Lu[t)v  I  B(v,v )  =  0,  t  >  0  ,  (3.12) 

»(0)  =  v0e  V  . 

Here  the  bilinear  operator  £?(•,■)  is  defined  using  the  Riesz  representation  theorem  as 

<  B(v,u),w  >v'xv=  6(w, «, t£?) ,  Vte  6  V  . 

The  evolution  Problem  4  can  be  derived  by  simply  applying  the  projection  operator  Ph 
onto  the  system  (3.1). 

Let  us  consider  the  following  integral  representation  for  the  solution  of  (3.12): 
Problem  5  Find  v  G  C([0,  oo);P)  such  that 

»(0  =  Z{t,0)v0-Jo  Z(t,r))B(v(r)),v{r]))dri  ,  (3.13) 

v(o)  =  vQ  e  v  , 

where  Z{ •,•)  is  the  evolution  operator  described  in  preceding  section. 

Theorem  3.3  Problem  5  resolves  Problem  4- 

Proof:  First  note  that  the  representation  (3.13)  in  Problem  5  formally  satisfies  the  evo¬ 
lution  form  (3.12)  in  Problem  4.  From  the  regularity  properties  of  the  evolution  operator 
Z{t,  0)  we  have, 

Z{t,0)v0  €  C([0,  oo);  P)  n  C1(0,  oo;  V') 

Hence  if  we  set, 

Y{t)  =  —  f  Z(t,rj)B(v(ri),v(r]))drj,  (3.14) 

Jo 

then  we  only  need  to  show  that  F(-)  £  C([0,oo);  V)  n  Cx(0,  oo;  V').  Let  us  first  show  that 
Y(t)  is  bounded  and  continuous  in  V. 


Note  first  that  B(-,  •)  £  C(V  x  V;  //).  In  fact  from  the  estimate  for  6(-,  •,  •)  in  (2.8)  with 
mi  =  2, m2  =  1 ,  m3  =  0  we  get 

|fc(u,v,tfl)|  <  co||«i|//i(n)||»||^{n)||to||^(n),  V«  £  V  ,  w  £  H  . 

Hence  by  the  Riesz  representation  theorem  we  can  write 

b(u,v,w)  =  (B(u,v),w)h*h  ,Vid  £  H  . 

By  setting  w  —  B[tt  t»)  £  //,  this  gives 

l|£(«,»)IU  <  coHMIHk- 


Thus  B(v,v)  £  C([0, 00); //)  for  v  £  C([0,  00):  V). 

Let  us  now  estimate  (3.14)  as, 

\\Y{t)\\v  ^  j0  llZ(^f?)ll£(H;V)||5(«(v),f(^))||//^  • 

Now  using  the  estimates  for  the  evolution  operator  Z{t,r)  obtained  in  the  previous  chapter 
we  get, 

ll^llc([O.T,);V)  <  Vl  (^2)  ll^llc(|0,Tj);V)  • 

Let  us  now  consider, 


Y(t  +  h)  ~Y{t)  =  -  [  \Z[t  +  h,r])  -  Z(t,r))}B(v(ri),v(ri))dri 

Jo 

rt  +  h 

-  Jt  Z(t  -f  h,rj)B(v{r]),v(ri))dr]  . 


Estimating  this  we  get, 

\\Y{t  +  h)~  y(0llv  <  { j*  II z(t  +  M)  -  z(t,i?)||£(W:ir)rfi? 

/«+* 

II Z(t  +  h,r))\\C{H-,v)dri}  || #(»,») ||c(o.r,;*)  • 


Again  using  the  results  of  the  previous  chapter  we  get 


||E(t  +  h)  -  y(OI|v  <  y2(<»||v||c(|0,T3);V) 
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with  y2 (t,h)  — *  0  as  h  — ►  0.  We  hence  conclude  that  F  G  C([0,  T2);F).  Now, 


Y{t  +  h)-Y{t)_  Z[t  +  h,t)-Z[t,t).  [* 


-( 


)  f  Z(t,rj)B(v(ri),v(Tj))dri 
Jo 


l  rt+h 

Z(t  +  M)^  Z(t,Tj)B(v(rj),v{T]))dT] . 


Taking  the  limit  h  — *  0  we  get  due  to  the  continuity  results  of  the  evolution  operator  Z( •,  •) 
obtained  last  chapter 


Yl(t)  =  -[eA  +  uAx  +  MOTO  -  B(v(t),v{t))  . 

Here  Y[  (•)  is  the  right  derivative. 

For  v  G  C([0,  r2);  V)  we  have  B(v,  v)  G  C([0, 7^);  H )  and  Y  G  C ([0,  T^);  V).  Since  A  is 
an  isomorphism  from  V  onto  V,  we  have  AY  G  C([0,  Tjj);  V')  and  \AX  +  Lu{t))Y  G 
C([0,  T^);  /f)  due  to  the  estimates  on  the  trilinear  form  &(•,•,•).  Hence  we  conclude  that 
the  right  derivative  exists  and  Y_[_(-)  G  C(0,T2;V).  From  this  as  before  we  conclude  that 
the  strong  derivative  F'(-)  G  C(0,T2\V').  Thus 

Y'(t)  +  (.AY  (t)  +  v>A{Y  ( t )  +  Z/y  (t)F(t)  =  0, 

which  implies  that  F (t)  is  a  solution  of  Problem  4  with  F  (0)  =  0.  This  proves  the  Theorem. 

□ 

The  existence  and  uniqueness  aspects  of  the  solution  can  be  established  using  the  meth¬ 
ods  used  in  [9].  Here  one  shows  that  for  a  fixed  time  interval  there  exists  a  neighborhood 
such  that  there  is  a  unique  solution  for  each  initial  data  in  this  neighborhood.  The  time  in¬ 
terval  can  be  taken  to  infinity  by  choosing  this  neighborhood  sufficiently  small.  (However, 
see  remark.) 

We  will  now  establish  an  important  result  regarding  the  dependence  of  the  solution  on 
the  initial  data.  Let  us  rewrite  Problem  5  in  the  form  J(w,r0)  =  0.  Here  the  map 

7(-,-):C(0,r2;F)xF->C7(0,r2;F) 
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is  defined  by 


7{v,v0)  =  v  -  Z(t,  0)t>0  +  M(v,v)  ,  (3.15) 

with  M(v,v)  =  [  Z{t,r])B(v(r]),v(T]))dT]  . 

Jo 

Here  the  bilinear  operator  M(-,-)  :  C(0,  T2\  V)®2  — +  C(0,  T2\V)  satisfies 

M(0,w)  =  M(w,0)  =  0,  VweC(0,T2;V ) 
and  Z(t,  0)  €  C(V,V)  for  fixed  t. 

The  following  theorem  is  a  consequence  of  the  analytic  version  of  the  implicit  function 
theorem  [5].  Results  of  this  type  for  the  conventional  Navier-Stokes  system  is  given  in 
Sritharan  [20], 

Theorem  3.4  For  fixed  T"  >  0  there  exists  a  neighborhood  of  the  origin  B2s  C  V  such  that 
for  v0  £  B2s,  there  exists  a  unique  solution  v  £  C(0,  T*\V)  to  the  problem  5.  Moreover, 
the  dependence  of  v  in  the  initial  data  Vo  is  Frechet-analytic. 

□ 


This  means  there  exists  a  map 

W(t,0--)  :B26CV  —>  Bls  C  C(  0,r*;T) 

such  that  v(t)  =  ty(t,0;t?0)  and  VV (t,  0;  •)  is  F-analytic  in  this  neighborhood.  That  is  r(t) 
can  be  written  as  a  power  series  in  the  initial  data  in  the  following  way: 


The  n-linear  maps 


v(t)  -=  W(t,0;vo)  =  Y,  ,v0;t). 

n>  1 

l/n(-  •  ■)  :  Bf"  — *  Big  continously  . 


(3.16) 


This  series  representation  converges  in  the  neighborhood  defined  above.  One  can  verify 
easily  that  )/i(i>o;0  =  Z{t, 0)vq.  We  finally  note  that  due  to  the  uniqueness  theorem, 


W(£,0;t?o)  =  W(Mi;^(*i,0; !>„)),  for  0  <  <!  <  t  <  oo, 
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and  VF(0,0;)  =  I. 

If  the  basic  flow  field  U(t)  is  T-periodic  then 

W (2T,  0;  t?0)  =  W{T,0;W{T,0;v0)) 

and  in  general 

W{nT,0\-)  =  W{T,0\-)n,  for  n>  1.  (3.17) 

As  noted  earlier  in  the  previous  section,  the  Frechet  derivative  of  the  solution  map  (DIF)(t, 
0;  •)  =  Z(t,0 )  satisfies  a  similar  relationship  as  in  (3.17). 

Remark:  The  nonlinear  semigroup  characterized  above  is  defind  (when  T'  =  oo)  only  in 
a  neighborhood  of  the  origin.  However  in  [15]  we  have  proved  the  global  existence  and 
uniqueness  of  strong  solution  in  V  using  other  methods. 

3.4  Local  Invariant  Manifold  Theorem 


Let  the  spectrum  of  the  mo nodromy  operator  Z{T,  0)  €  C{V',V)  splits  into  two  disjoint 
sets  Ofj  and  os  such  that  a(Z(T,  0))  =  ov  U  as  and 


bs  =  sup  |A|  < 

Ae<7S 


inf  | A |  =  b 

Aec j„ 


-l 

u  • 


Let  P{j  and  Ps  be  the  spectral  projectors  defined  by  the  Dunford’s  integrals, 


(3.18) 


Pv  =  ~[  R(\-,Z(T,0))d\ 

2 7ri  Jr u 

and  Ps  =  —  f  R{X-Z{T,0))dX. 

27u  J rs 

Here  R(X;Z(T,  0))  is  the  resolvent  operator  and  ry,r5  encircle  ou,os  respectively.  Note 
that  Ps  +  Pu  =  1  ,  PsPu  =  PuPs  and  Ps,Pu  commute  with  Z(T,  0). 


Theorem  3.5  (The  invariant  cone  theorem)  Let  the  basic  flow  be  T-periodic  in  time 
so  that  the  solution  map  satisfies  (8.17). 

(i)  If  the  spectrum  of  the  monodromy  operator  Z(T,  0)  lies  inside  the  unit  disc  (spectral 
radius  <  1 )  then  the  basic  periodic  solution  is  (locally)  exponentially  stable:  there  exists 
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p  >  0  such  that  V  v0  G  5^,(0)  C  V,  W(t,0;w0)  — >  0  exponentially  in  the  norm  of  V . 

(ii)  Let  the  spectrum  of  Z(T,  0)  satisfy  (S.18)  with  by  <  1,  then  there  exists  a  double  cone 
\  c  V  and  a  ball  Bf( 0)  C  V  such  that  Vt>0  G  B^(0)  n  X'\{0},  there  exists  n  £  H  for  which 
\V(nT,  0;  r0)  J[ v  >  <5-  That  is  the  basic  solution  is  Lyapunov  unstable.  Here  double  cone  is 
defined  by 

K  ~  {u  G  V  such  that  j|Ps®||v  <  -)f[[P(/t>||v,  7  >  0}.  (3.19) 


□ 

Theorem  3.6  (The  invariant  manifold  theorem)  Let  bs,bu  <  1.  Then  in  a  neigh¬ 
borhood  Br(0)  C  V ,  there  exists  two  unique ,  analytic  manifolds  Ms  and  Mu  which  are 
respectively  the  graphs  of  the  maps,  <j>s  :  PsV  — *■  PuV  and  fiu  ■  PuV  PsV  ■  maps 

4>s  and  4>u  are  analytic  with, 

1.  cf>u(0)=4>s{0)  =  0; 

2.  D(j)u{ 0)  =  D(j>s{ 0)  =  0  (tangency  condition;) 

3.  manifolds  Ms  and  Mu  are  locally  invariant  under  the  solution  map  \V(T,  0;  ■) 

W{T,  0;  Mu  n  BT  (0) )  C  Mv  and  W{T,  0;  Ms  n  BT  (0) )  C  Ms ; 

4.  stable  manifold  Ms  satisfies 

Ms  H  Br{ 0)  =  [v  G  Br(0)  such  that  V  n  >  0,  \V{nT,  0;  v)  G  Br{ 0) 

and  — >  0  as  n  — ►  00}; 

5.  Unstable  manifold  Mv  satisfies, 

Mv  n  Br(0)  ={»€  Bt( 0)  such  that  VV(T,  0;  • )nv  is  defined  Vn  <  0 
and  tends  to  zero  as  n  — +  ~oo}; 

6.  if  v  Ms  then  there  exists  6  >  0  and  p  G  /V  such  that,  \\\V(pT,  0;  v)||v  >  6; 
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7.  dist(M[/,  W (T,  0;  t>))  <  dist(A fu,v)  for  v  G  i?r(0)  (exponential  attractive  property 
of  the  unstable  manifold;) 

8.  dist(Ms,  W (T,  0;  t?))  >  dist(A/5,v)  for  v  G  BT{ 0)  (repelling  property  of  the  stable 
manifold.) 

Proofs  of  Theorem  3.5  and  3.6  are  similar  to  those  for  the  conventional  Navier-Stokes 
system  and  can  be  found  in  detail  in  [20]. 

□ 
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Chapter  4 

Global  Invariant  Varieties  (Inertial 
Varieties) 


Foias,  Sell  and  Temam  proposed  in  [6]  the  concept  of  inertial  varieties  for  certain 
class  of  semilinear  evolution  equations.  In  this  section  we  will  study  the  existence  of  such 
global  invariant  varieties  modelled  on  the  invariant  subspaces  of  A.  These  manifolds  will 
be  invariant  to  the  action  of  the  W[t,  0,  •).  In  section  4.5  we  will  extend  the  general  theory 
in  [6]  to  analyze  the  regularity  of  the  inertial  manifolds.  Let  us  first  obtain  certain  overall 
bounds  for  the  solution  in  various  norms  and  show  that  the  dynamics  is  characterized 
by  a  compact  global  attractor.  The  global  manifolds  to  be  constructed  will  contain  this 
attractor. 

4.1  Overall  Bounds  for  the  Solutions 

We  can  get  a  weak  formulation  from  the  system  (2.1)-(2.4)  by  taking  duality  pairing 
with  to  6  V, 

du 

<  — ,  to  >  +t(Au,  Ato)  +  i/(Vt»,Vto)  +  6(u,  u,  to)  =<  /,  to  >  V  to  G  V  (4.1) 

C *  V 

We  will  first  consider  the  energy  estimate  by  setting  to  =  u  and  using  the  fact  that 
b(u,  u,  a)  =  0, 

\jt M2  +  eM2  +  v\ V“|2  =</,“>•  (4.2) 
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Using  the  Poincare’s  Lemma  |Vti|2  >  r/Jul2  and  the  fact  ||tt||2  >  Pi\u\2  we  get 

^-\u\2  +  ( £ pi  +  2urjl)\u\2  <  Miiui. 
cLt  € 

Here  we  denote  rj j  as  the  smallest  eigenvalue  of  Stokes  operator  A\.  We  then  set  a  = 

ll/ll2' 

epi  +  2 vr)i  and  pi  =  -  to  obtain 

cot 

|u(i)|2  <  |u(0)|Vaf  +  p2q  [  1  -  e~at\,  a  >  0.  (4.3) 

We  note  that  in  the  case  of  |«(x,  0)  |  >  p0 ,  the  energy  estimate  in  (4.3)  show  that  |u(<)  |  is 
bounded  by  a  monotone  decreasing  exponential  funtion.  On  the  other  hand,  if  |u(x,  0)|  < 
Po,  then  |«(x,t)|  <  po  for  Wt  €  R+ ■  Hence  for  any  ball  Br  —  («(0)  E  H\  |«(0)|  <  R},  there 
is  a  ball  Bro  in  H  centered  at  origin  with  radius  r0  >  pQ  such  that 

W  {t,  0;  Br)  C  Bro,  for  t  >  t0{BR)  =  -In  ^4-  (4.4) 

a  r0  Po 

The  ball  Brg  is  said  to  be  exponentially  absorbing  and  invariant  [10,11,6]  under  the 
action  of  the  map  W ( t ,  0;  •). 

Let  us  now  proceed  to  get  other  estimates.  Notice  that  from  energy  estimate  if  u0  6  BR 
and  t  >  t0(BR),  by  integrating  (4.2)  from  t  to  t  +  1, 

/t+i  1  1 

\\u\\2dr<^[r20+i7\\f\\2vl)  =  ao.  (4.5) 

To  prove  uniform  bounds  on  different  norms  we  use  the  uniform  Gronwall  inequality  [6,15]: 

Lemma  4.1  (Uniform  Gronwall  Inequality)  Let  g,h,y  be  three  positive  locally  inte- 
grable  functions  for  t*  <  t  <  +oo  which  satisfy 

—  <  gy  +  h  V  t  >  t*  and 

dt 

/£+l  rt+ 1  rt+l 

g(s)ds  <  cti  ,  J  h(s)ds  <  <*2  ,  j  y{s)ds  <  0:3  , 
for  all  t  >  t* ,  where  ai,a2,a3  are  positive  constants.  Then 

y{t  +  1)  <  («2  +  a3)exp(al)  ,  V  t  >  t\  (4.6) 
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Let  us  now  suppose  that  /  <E  H  and  set  w  =  Au  in  (4.1), 

+  e\Au\*  +  Mv^2«|2  <  Ci|u|1/2||«||  |Au|3/2  +  |/|  |A«|. 

1  at 

Here  we  have  used  the  estimate  in  (2.10)  for  trilinear  form  b(u,u,  Au).  Using  the  Young’s 
inequality,  we  obtain  the  following  differential  inequality  after  dropping  the  positive  term 
i/|VA2u|2 

^•||u||2  +  e\Au\2  <  c;|«|2||«||4  +  (4.7) 

We  showed  earlier  that  any  solution  will  enter  an  absorbing  ball  Bru  in  H  for  t  >  t0(BR). 
Thus 

This  result  together  with  the  estimate  (4.5)  allows  us  to  apply  the  uniform  Gronwall 

2I/|2 

inequality  with  y  =  ||«||2,g  =  c'1(r2)||u||2  and  h  =  - .  we  get 

2l  /I2 

HOII2  <  (^0  H - exp(c'i  r2  o0),  for  t  >  i0{BR)  +  1, 

where  t0(BR)  is  given  in  (4.4).  This  means  that  there  exists  an  absorbing  ball  Bri  of 
radius  rq  in  V  such  that  all  the  solutions  will  enter  this  ball  after  certain  time: 

W(t,0\BR)  c  Bri,  for  t  >  t0(BR)  +  1. 

Since  Bri  is  compact  in  H  we  conclude  that  W ( t ,  0;  •)  maps  bounded  sets  in  H  into  compact 
sets.  Note  that  it  follows  from  (4.7)  that  for  t  >  t0(BR)  +  1, 

rt+l  t  2lfi2 

/  |A«|2dr  <  -  (ci  r2r4  H - ]  =  cq.  (4.8) 

Jt  e  e 

Now,  let  us  set  w  ~  A2u  in  (4.1)  to  get 

\jt\Au\2  +  cllAult2  -  u\Au\ llAull  +  CziMl  Hwl!  IIA“II +  ll/ll  WAuW- 

Since  any  solution  will  be  absorbed  by  balls  Bro  and  Bri  after  time  t  >  t0(BR)  +  1,  this 
becomes 

d  |  .rt  l2z^2|  ,2  42c2/2w2\  12 1|  _  ||  2 

-|A«|2  <  - |Au|2  +  — ^(r2)(r?)  +  -||/||2. 

at  e  e  e 
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Hence,  we  can  apply  the  uniform  Gronwall  inequality  again  to  conclude  that 

12  10  f/2 

|A«(t)|2  <  [<7i  +  —  {c'2r20rl  +  ||/]|2)  ]  exp(— — ),  for  t  >  t0(BR)  +  2. 

That  is  ,  there  exists  an  absorbing  ball  Br 2  of  radius  r2  in  D(A )  such  that  all  the  solutions 
will  enter  this  ball  after  time  t  >  t0(BR)  +  2: 

W(t,0]  BR)  C  Brj,  for  t  >  t0(BR)  +  2. 

Note  that 

closure  of  W{t,0;BR ) 

t>to(BR)  +  l 

is  compact  in  H.  Now,  the  compactness  of  the  operators  VT(t,0;  •)  in  H  implies  that  there 
exists  a  compact  attractor  A  which  attracts  every  bounded  sets  in  H.  In  fact,  A  is  the 
global  attractor  for  the  operators  W  ( t ,  0;  ■)  and  it  is  also  the  cu-limit  set  of  absorbing  set 
Br2,  i.e.  A  =  u(BT2).  This  means  if  we  denote  W(t,  0;  Br2)  =  Br2(t )  then 

A-n4u*.M). 

T> 0  \t>T  J 

Note  that  the  global  attractor  A  must  be  contained  in  the  absorbing  balls  in  H,  V  and 
D{A): 

icj,nsrir,j,. 

In  addition  one  can  show  that  if  W(t,0\-)  is  injective  then  in  A,  iy(£,0;-)  will  be  defined 
for  all  t  6  R  [19].  Such  a  result  for  2-D  conventional  Navier-Stokes  equations  has  been 
proven  by  Ladyzhenskaya  [10].  In  the  appendix  we  prove  the  time  analyticity  of  the  map 
IT(t,0;  ■)  which  implies  injectivity. 

Remark:  Note  that  if  we  assume  /  €  H,  then  from  the  energy  estimate  we  also  get 

4lul2  +  (2e^i  +  <  —  • 

dt  vr/x 

I/I2 

By  setting  a'  =  2e^i  +  urjx  and  p’S  =  - 1 —  we  obtain 

urjia' 

]«(<)|2  <  |u(0)|2e-a,t  +  p'02\l-  e~a,t  ],  a' >  0.  (4.9) 
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Notice  that  the  above  estimate  is  uniform  in  e.  Let  us  denote  W{t,0;  Br>a)  =  Br>g(t)  and 
note  that  (4.9)  implies  that  the  absorbing  ball  Br'o(t)  will  shrink  as  t  increase.  That  is 

Br>o(t2)  C  Br<g{t 0,  for  t2  >  h. 

Note  that  from  energy  estimate  we  can  also  deduce, 

rt+ 1  1-1 

/  |Vu|2dr  <  -  [tq  H - |/I2  ]  =  cto  for  t  >  tQ(BR<). 

Jt  u  urji 

Let  us  now  set  w  —  Axu  in  (4.1), 

i-^-|Vu|2  +  e|VA2u|2  +  iy\Aiu\2  <  ci|u|1/2|Vu|  |>litx|3'/2  +  |/|  (4.10) 

2  at 

In  the  above  we  have  used  the  following  estimate  for  trilinear  form  b[u,u,  Axu)  which  is 
valid  only  for  n  =  2: 

|6(U,  U,  Axil)  <  C!  |t*|1/2|  Vtx||Ai«|3/2. 

Using  the  Young’s  inequality  in  (4.10)  and  then  dropping  positive  terms  e|Vj42u|2  and 
u\Aiu\2,  we  get 

—  |Vu|2  <  c'j|«|2IVu|4  + 
dt]  1  u 

By  applying  the  uniform  Gronwall  inequality  again  to  obtain 

|Vu(<)|2  <  (o0  +  exp(c'x  r'02  o0),  for  t  >  t0(BR .)  +  1. 

Note  that  the  estimate  is  again  uniform  in  e.  This  means  that  there  exists  an  absorbing 
ball  Br't  in  V*  =  {u  6  Hq(U)]  divu  =  0}  such  that 

W(t,0;  BR')  C  for  t  >  t0(BR< )  +  1. 

Let  us  denote  the  global  attractor  At  as  w-limit  set  of  absorbing  set  Br i .  This  means  if 
we  denote  VU(t,0;Br')  =  B^ft)  then 

A<  =  fl  cl\  U  Br\  (r) 

r>0  \t>T 
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Note  that  since  W(t,  0;  •)  maps  a  bounded  set  in  H  into  compact  set  in  H,  At  is  compact 
in  II.  (This  is  due  to  the  compactness  of  embedding  V *  in  H ).  All  the  estimates  used 
above  are  uniform  in  e.  Such  results  will  be  useful  in  establishing  the  limit  of  Af  to  the 
attractor  A*  of  the  conventional  Navier-Stokes  equations  as  e  — *  0.  Special  cases  of  such 
limit  solutions  were  established  in  [15]. 


4.2  Formulation  of  Inertial  Varieties 


Let  us  now  consider  the  system  (2.l)-(2.4)  as  an  equation  of  evolution  in  Hilbert  space 


V : 


du 

dt 


+  eAu  +  vAiU  +  B(u,u)  =  f,  t  >  0, 


(4.11) 


u(0)  =  u0  6  V. 

In  the  sequel,  we  will  assume  /  G  H . 

In  the  previous  section  we  observed  that  although  absorbing  balls  exists  in  II ,  V  and 
D{A)  spaces,  the  positive  invariance  property  can  be  established  only  for  the  Lf-ball  Bru. 
However,  we  are  interested  only  in  the  dynamics  in  a  neighborhood  of  the  attractor.  We 
will  thus  devise  a  method  to  restrict  our  investigation  to  the  dynamics  inside  the  absorbing 
ball  Bri.  We  will  use  a  smooth  cut-off  funtion  0  :  R+  — ►  [0,1]  and  #r,(r )  =  0(r/ri)  such 
that 

f  0(£)  =  1,  for  0  <  £  <  1 
6{£)  =  0,  for  £  >  2 
l  TOI<2,  for  £>0, 

to  modify  the  equation  (4.11): 


(4.12) 


where  R{u)  =  i/AiU  +  B(u,u)  -  /. 

The  absorbing  property  of  this  modified  equation  can  be  shown  by  taking  inner  product 
with  Au  to  (4.12).  Let  us  consider  the  solution  starting  outside  the  ball  B2ri  in  V,  i.e. 
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||u||  >  2 jq.  Then  we  obtain  the  following  simple  differential  inequality  since  0ri(||u||)  =  0 

^ll«H2  +  (fMi  +  *"?i)IM|2  <  “||«||2  +  el^u!2  +  HV^2«|2  =  0. 

Here  we  have  used  the  fact  that  /^Uull2  <  |Au|2  and  t7j||u||2  <  |VA2tt|2.  This  leads  to 
||«(0I12  <  ll“(0)||2e~5£,  with  (3  =  2(c/£t  +  urji). 

Notice  that  if  ||ti(0)||  >  r3  with  r3  >  2ri,  then  the  solution  orbit  will  converge  expo¬ 
nentially  to  a  ball  Bri  of  radius  r3.  If  |]tt(0)j|  <  r3,  then  solution  orbit  will  stay  inside  the 
ball  BT3  and  never  leave  this  ball.  That  is,  there  exists  an  exponentially  absorbing  and 
invariant  ball  Br3  in  V  for  this  modified  equation  (4.12).  In  addition,  from  the  definition 
of  the  cut-off  funtion  6,  we  have  0ri(||u||)  =  1  for  ||u||  <  rq  (i.e.  «  £  Bri).  This  means  that 
the  original  equation  (4.11)  and  the  modified  equation  (4.12)  are  identical  in  the  neigh¬ 
borhood  of  the  global  attractor.  Thus  tm  the  dynamics  of  (4.11)  are  exactly  represented 
by  system  (4.12)  after  a  certain  time. 

We  shall  define  the  solution  operator  Uo  — ►  t i(t)  associated  to  the  modified  equation 
(4.12)  as  Wm(t,0\-).  The  operator  Wm(t,  0;  •)  is  defined  for  all  t  >  0  and  all  u0  £  H. 
The  existence  of  the  absorbing  invariant  ball  BT3  C  V  implies  that  Wm(t,0;u0)  £  BT3  for 
t  >  Tc(u0)  for  any  initial  condition  u0  £  V.  We  will  thus  construct  the  inertial  variety 
At  with  a  compact  support  in  Br3.  According  to  Foias,  Sell  and  Temam  [6],  the  formal 
definition  of  inertial  variety  is  as  follows:  A  set  X  C  V  is  called  an  inertial  variety  of 
(4.12)  if 

1.  At  is  a  finite  dimensional  Lipschitz  manifold; 

2.  At  has  a  compact  support  and  is  invariant  under  the  solution  operator  Wm(t,0;-), 
i.e.,  Wm(t,0;  W)CX  for  all  t  >  0; 

u.  all  orbilo  of  the  solution  of  equation  (4.12)  in  V  are  attracted  exponentially  to  At. 


Remarks: 


(a)  We  begin  with  this  formal  definition  given  in  [6]  and  then  study  the  smoothness  of 
such  manifolds  in  section  4.5. 

(b)  Note  that  the  condition  3  implies  that  Ac  X.  In  fact,  if  uA  G  A,  then  since  A  is 
its  own  w-limit  set,  for  each  t  >  0,  3tiA  G:  A  such  that  uA  =  Wm(t,0;vA).  Now 

dist(«A,  X)  =  dist(Wm(t,0;wA),  >1)) 

<  e~atdist(t>A,  X), 

where  a  is  a  positive  constant.  Taking  t  — >  oo,  we  conclude  that  dist(uA,  X)  has  to  be 
zero.  Since  X  fl  Bri  is  closed,  we  conclude  that  uA  G  X. 

Let  us  now  define  p  =  P^u  and  q  =  QNu  for  u  £V.  By  applying  the  projection  P n  and 
Qn  to  the  equation  (4.12)  we  obtain  evolution  equations  for  P^V  and  Q^V  respectively. 

+  eAp  +  PnF[u)  =  0,  (4.13) 

at 

^  +  eAq  +  QNF(u)=0,  (4.14) 

where 

f(«)  =»r.(ll«ll)B(«)- 

Let  us  construct  X  as  the  graph  of  the  Lipschitz  function  $  :  P^V  — ►  Q^V  with  $ 
belonging  to  the  class  Hbt  defined  below. 

Definition  4.1  Let  Hbl  be  the  space  of  Lipschitz  maps  <J>  from  PNV  into  Q^V  which 
satisfy 

(i)  ||*(p)||<&,  b>0,  V  p  G  PNV, 

(ii)  ||$(p,)  -  $(p2)||  <  /  ||px  -  p2||,  Vp,,p2  G  PNV, 

(iii)  supp  $  C  {p  G  PNV  :  ||p||  <  r3}. 

Note  that  LIbl  is  complete  with  the  metric 

dist  (*,*')  =  11$  -*'|k,=  sup  ||*(p)-*'(p)ll- 

vePNv 
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Let  $  E  Hu  and  p0  E  PjvV,  then  p(t )  will  be  determined  by  the  following  initial  value 


problem: 


+  eAp  +  PNF(p  +  $(p))  =  0 
at 

P  =  P(t)  =  p(t;*,p0) 


(4.15) 


l  P(  0)  =  P0 

Let  us  first  note  that  the  modified  nonlinear  term  F(u)  is  globally  bounded  and  Lips- 


chitz. 


Lemma  4.2  For  $  6  Hb,i  and  pt,p2  G  PnV ,  we  have 


(i)  \F(u)\<dl, 

(ii)  |F(«1)-i;,(«2)|<^(l  +  /)||Pl-P2l|, 
where  d2  —  2 rx1dl  -f  v  +  4c0r1  and  dx  =  2 urx  +  4c0r?  +  |/|. 

Proof:  (i):  Note  that  we  have  the  following  estimates  for  A i  and  J3(u,r) 

|4i«|  <  ||«||,  VuGL 
|J3(«,w)|  <  c0||«||  ||r||,  Vu,rGF. 


(4.16) 

(4.17) 


(4.18) 

(4.19) 


This  gives 


\F{u)\ <  0ri(!!«li)  [v\Axu\  +  \B(u,  u)|  +  |/| 

<  MINI)  [HMI  +  colMI1  + 1/|] 

<  2 i/rx  +  4corj  +  |/|  =  dx. 


(4.20) 


where  dx  is  independent  of  time, 
(ii)  We  write, 


R(ux)  -  R{u2)  =  ^Ai(t»i  -  u2)  +  B(ui,ui  -  u2)  4-  5(ut  -  u2,tt2)- 


Using  (4.18)  and  (4.19)  we  get 


\R{ux)  -  iZ(t»2)|  <  HIU1  -  «2||  +  C0||Ui||  ||«!  -  U2||  +  Co||«l  -  «2  II  l|«2| 


33 


That  is, 


l-R(ui)  -  R[v*)\  <  {v  +  co||«i||  +  c0||«2||)||«i  -  «2||.  (4.21) 

Let  us  investigate 

F(u,)  -  F(u,)  =  *ri(ll«.ll)K(«i)  -  #ri(||«,||)K(«1) 
for  the  following  three  cases. 

1-  ||«iii>IM  >  2rt, 

2 .  1 1 XJL j  1 1  2Tj  ^  1 1 1^2 1 1  ||  ^2  ||  2rj  ^  IM. 

3.  ||ui||,  ||u2||  <  2rx. 

In  the  first  case,  since  both  ^ri(||tii||)  =  0  and  ^ri ( || ti2 1| )  =  0  we  have 

F(«i)  —  F(u2)  =  0.  Now  if  Huxll  <  2rx  <  ||«2||  then  since  0ri(|jt/2||)  =  0  we  have  ^(ui)  — 
F(u2)  =^1(||tt1||)i2(ui).  Thus 

<  |^ri(||«i||)  -  tfri(i|t»2||)|  |J2(t*i)| 

<  2rf1|  Huill  -  ||u2||  |  |i2(«1)| 

<  2rI_1di||«i  -  u2||,  (4.22) 

where  di  =  2 ur\  +  4corJ  +  |/|.  A  similar  result  can  be  obtained  for  ||u2||  <  2 rx  <  ||«i||.  For 
the  last  case  we  have 

H «.)  -  F(«l)  =  [  *-,(ll«.||)-«r.(ll«>ll)  ]«(«.)  +«r.(ll«»||)i  *(«.)-  *(«»)]• 

From  (4.21)  and  (4.22),  we  get 

|F(«i)  -  F(«2)|  <  {  2r1'1d1  +  v  +  4c0ri  }  ||ui  -  u2|| 

<d2||t*i-«2||,  (4.23) 

with  d2  =  2rf  1d1  +  u  +  4c0ri. 
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Let  us  estimate 


«i  -  «2  -  Pl  -  p2  +  [$(Pl)  -  $(p2)]. 

as, 

ll«.  -  «nll  <  (1  +  OIIp,  -  p,ll- 

Combining  this  with  (4.23)  gives 


l-F(ui)  -  F(u2)|  <  d2  (1  +  l)\\pl  -  p2|[. 


□ 

In  consequence  of  (4.17),  we  have 

\P"F(Pi  +  *(Pl))  -  PNF(p2  +  $(p2))|  <  d2  (1  +  /)||Pl  -  P2||,  Pl,P2  €  PNV. 

This  means  P„F  is  a  mapping  from  PNV  into  PNH  and  satisfies  global  Lipschitz  condition. 
From  a  uniqueness  theorem  for  evolution  equations  in  finite  dimensions  [2,3],  there  exists 
at  most  one  solution  p(t)  =  p(t;$,p0 )  for  the  Cauchy  problem  (4.15). 

Let  us  now  integrate  (4.14)  from  £  to  t  to  get 

fl(t)  =  QNS(t  —  £)?(£)  +  QNS(t  —  t)QnF{u{t))(It. 

Notice  that  (4.16)  implies  that  QNF(p+  $(P))  is  a  bounded  operator  such  that  QNF  : 
PnV  ->  QnH.  Hence,  there  exists  a  unique  solution  q{t)  =  ?(<;$,  Po)  that  stays  bounded 
as  £  ~ *  -oo.  This  can  be  seen  as  follows:  from  Lemma  3.3 

||Qw5(t  -  f)?U)ll  <  ->  0  as.  £  ->  -oo. 

Also  it  follows  from  (4.16)  that 

I"'  ||C?JV5(t-r)QA,F(u(7-))||dr 

*  —  OO 

—  j  ^  II <5 at-S (i  -  ’')!!£(Qiv//;QjvVq  \QnF(u(t)) \dr 
—  II  QNS(t  -  T)\\£(QNff.QNV)dT 

<  dxrlnNll\  for  t  e  R.  (4.24) 
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We  now  take  £  =  —  oo 


q{t)  =  -  [  QNS{t  -  t)QnF[u{t))(1t.  (4.25) 

J  -oo 

Here  S(£)  is  the  holomorphic  semigroup  generated  by  —eA  described  in  section  3.1.  Since 
q(t)  is  continuous  and  bounded  for  all  t  G  R,  we  choose  t  =  0  so  that  p0  G  P^V  will  be 
related  to  ^(0)  =  g(0;$,p0)  by 

?(°)  =  -  [  QnS(-t)QnF(u(t))(It. 

J  —  OO 

We  thus  seek  the  (unique)  uniformly  bounded  solution  of  (4.11)  on  (— oo,0]  with  u  = 
p  +  $(p).  We  will  now  define  a  mapping  T  :  p0  — >  g(0)  as 

(r$](p0)  :=  -  f  Q p/S(  —  t)Q n F(u(r))d,T,  G  Hbl,  (4.26) 

J  -oo 

where  u(r)  =  p(r;<I>,p0)  +  $(p(r;  4>,p0)).  Notice  that  the  condition  that  u  =  p0  +  g(0) 
belonging  to  At  is  equivalent  to  the  existence  of  a  fixed  point  for  the  map  T : 

9(0)  =  *(p0)  =  r$(p0),  V  p0  G  PNV. 

4.3  Inertial  Manifold  Theorem 

We  will  establish  the  existence  of  inertial  manifold  using  contraction  mapping  theorem. 
For  this  purpose  we  must  prove  the  following: 

Lemma  4.3  For  $  G  Hb  l,  if  pjj2+l  —  >  e~ld3  and  l  <  2  then  the  T  maps  Hb  i  into 

itself: 

T  :  HKl  — »  Hbth 

where  d3  =  2d2(l  +  0^  1  an< 1  ^  *s  the  Lipschitz  constant  for  $ 

Proof:  We  first  note  that, 

Lemma  4.4  For  $  G  Hb  i,  we  have 

supp C  (p  G  P^V  :  ||p||  <  r3},  where  r3  >  2 rt. 
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Proof:  We  need  to  show  for  ||p0||  >  r3  we  have  T<J>(p0)  =  0  for  all  $  £  lh,i-  Notice  that 
u  =  Psp  +  Qjv$(p)  since  Pf,  =  Pff  and  Q2n  =  Qn-  Now  A'^u  ~  PNAx/2p  +  Q N Ax!2${p) 
since  P^,QN  commute  with  A1/2.  Thus,  ||u||2  =  ]Jpj|2  +  ||<F(p)||2.  Hence  for  J]p||  >  2rx  we 
have  ||u||  >  ||pjj  >  2 rt.  This  implies  0ri(||u||)  =  0  for  ||p||  >  2ry. 

Let  us  consider  the  initial  point  outside  the  ball  Brz  such  that  ||p 0[(  >  r3  >  2r\.  We 
can  then  find  a  time  t  such  that  ||p(f)||  >  2r1}  for  t  <  t  <  0  and  hence  0ri  ( j| u 1 1 )  =  0.  Thus 
(4.15)  becomes, 


%  +  <A,  =  0 

dt 


l  P(0)  =  Po- 

By  taking  inner  product  with  Ap  for  above  system  and  using  the  fact  \Ap\2  >  Mi l| p|j 3 


we  get 


^llpll2  A  <  0. 


This  gives  us 


I|P(0)||  <  \\p{r)\\e^T  <  ))p(r)l  r  <  0. 


Since  j|p0||  >  r3  >  2 rx,  we  have  H p(t) ||  >  r3  >  2ri.  This  means  that  u(r)  always  stays 
outside  the  ball  BT3  for  all  r  <  0.  Therefore,  0r,  (])«]))  =  0  for  all  r  <  0.  This  gives  us 
F(u)  =  0  and  thus  T<b(p0)  =  0. 


□ 


Let  us  show  that  T$  is  globally  bounded  in  V . 

Lemma  4.5  If  p0  £  P^V ,  then  (T$j(p0)  £  QnV  and 

l|r*(p0)ll<&, 

with  b  -----  2d\t~le~xl2 Pxl\ . 

Proof:  From  the  definition  of  the  map  T  it  is  obvious  that  T$(p0)  £  Q^V  and 

l|r*(p0)||<  f  ||Qyv5(-r)||c(0;vK;QNv)l<3NF(«)|dr. 

J  -oo 
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Recalling  the  estimate  for  from  Lemma  3.3  (with  a  =  1/2)  we 


[  \\QnS(-t)\\z{qnH.qnV)<It  <  2e  xe  ^Vw+i- 

J  —  oo 


It  follows  from  (4.27)  and  (4.16)  that 


(4.27) 


|r$(p0)||  <  2d1e-1e-1/V)v1+/12  =  b. 


(4.28) 


Let  us  now  prove  the  global  Lipschitz  continuity  of  the  map  [T4>](-). 


Lemma  4.6  If 


1  /2 

tiff  —  e(fJ.N+ i  —  P-n)  ~  d2(l  +  OMn  >  0, 


then  for  $  £  Hb  l  and  for  p01,p02  £  toe  have 


||  [r$](Pol)-[r$](p02)  ||  </'||Pbi-Pb2ll, 

with  /'  =  d2(l  +  l)n~Nll\  [e_1  +  (e  -  cr^^)_1]  e~1/2exp 

and  oN  =  ■  ^'V-,  Cn  =  e  +  d2(l  +  l)p,^2 . 
Pn+  i 


2e 


(4.29) 


(4.30) 


Proof:  Let  $  be  a  fixed  element  in  f/jj.  Let  Pj  =  Pj(i)  and  p2  =  p2(t)  be  two  different 
solutions  of  the  initial  value  problem  (4.15)  with  Pj(0)  =  Poi>P2(0)  =  P02  respectively. 


+  eApy  +  P^P(ui)  =  0 


(4.31) 


Pi(°)  =  P  01 


~~  +  (A  p2  +  PatP(u2)  —  0 
at 


(4.32) 


l  P2(°)  =  Po2, 

where  Ui  =  Pj  +  $(px)  and  «2  =  p2  +  $(p2).  By  subtracting  (4.32)  from  (4.31),  we  get 


-jj  +  eAp  +  JW(«i)  -  PtfF(u2)  =  0, 
dt 


(4.33) 
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where  p(t)  =  p^t)  -  p2(t )  and  p(0)  =  p01  -  p02. 

VVe  take  inner  product  of  Ap  with  (4.33)  and  use  estimate  (4.17)  to  get 


p||2  +  e|Ap|2  <  |Pjv[F(u!)  -  F(u2)}\  \Ap\ 
<  d2(l  +  /)||p|| \Ap\. 


By  using  the  fact  |Ap|2  <  Miv || p||2  we  obtain 

IIpII^IIpII  >  -  [ePN  +  d2(l  +  l)pj2}  ||p||2 


which  implies 


^IIpII  +  PnZn  ||p||  >  0,  where  £jv  =  e  +  rf2(l  +  l)pN1/2. 


This  leads  to 


||Pi(0  -  P2(r)ll  <  ||Poi  -  Po2ilexP{-PN£jvT},  V r  <  0.  (4.34) 

Now,  using  the  Lipschitz  condition  on  F(u)  given  in  (4.17),  we  can  estimate 

lir*(Poi)  -  r«(pra)||<  f  ||<?ws'(-T)||£W„i,1,„p)|«w[F(«I)  -  f («,)]!* 

J  -oo 

<  d3{l  +  1)  I  ||Q^5(-r)||£((3wi/;(3jvV)||p1(7-)  -  P2(r)||dr 

J  —  OO 

By  applying  result  (4.34)  we  get 

i  /•>  f-l/^enN  +  i 

||r*(Poi)-r*(pte)||  <  d2(l  +  OIIPoi  -  P02II  Willi  /  exp [6NT\dT 

J  -CO 

/  1  x1/2  rO 

+  (5—)  /  |r|“1/2exp  [-nN^NT}dT}t 

\6£C/  »  — 1/2  (fiftf  +  i 


where  6^  =  e{pN+i  -  pN)  -  d2{  1  +  l)p 


1/2 
jv  ■ 


Note  that 


1/2  r-'/i^N+i 

Mn+i  / 

«/  —  OO 


e*wTdr  <  e  (e  ~  1 


°n(n 
e  2< 


) 


with  >  0  and  =  pN/pN+i.  Also 


(— )1/!  f 

\2eeJ  J-i/2€vn+i 


|r|  1/,2e  tiN^NTd.T  <  e  *e  1^2pAf1/12 


(4.35) 


(4.36) 
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We  thus  get, 


l|r#(Poi)  t<j>(p02)||  <  i'  |jpoi  ~  P02II) 

with 

l'  =  d2(  1  +  Z)Mw+?  [c_1  +  (c  -  &n £n ) _1  ]  e_1/2  exp  ~i~- 

□ 

Now  in  order  for  the  transform  T  to  map  Hbl  into  itself,"  we  must  have  V  <  l.  An 
elementary  calculation  shows  that  sufficient  conditions  for  l'  <  l  are 

e~ld3  <  Mwh,  with  d3  =  2d2(l  -f  /)/_1  and 
ea'J2  +  d3  p-Nxll  <  with  /  <  2. 

Combining  above  two  results, 

Hn+i  ~  pjy2  ^  e_1^3>  with  d3  =  2d2(l  +  l)l~ 1  and  /  <  2. 

□ 

Lemma  4.7  For  $  G  1/  >  e_1d4  and  /  <  2  then  the  transform  T  is  a  strict 

contraction  on  Hbi.  Here  dA  =  2d2(2e-1/2  +  /). 

Proof:  We  will  first  show  that  for  $1($2  £  an(^  Po  £  -PjyV", 

II  [r$i(p0)  -  r$2(p0)]  ||  <l||*1-«2||. 

Consider  two  elements  and  $2  in  Hb  i.  We  take  U\  =  Pi  +  $i(Pi)  and  tt2  =  p2  +  $2(p2) 
with  same  initial  condition  p0.  Then  analogous  to  (4.33)  we  get 

^  +  eAp  +  -  PnF(u2 )  =  0 

p(0)  =  0. 

By  taking  inner  product  with  Ap  and  again  using  the  fact  |Ap|2  <  p^||p||2  we  get 

4-IIpII  +  Pw6v||pII  >  -diPN2  ||$i  -  $:|| 

,  at 

.  p(0)  =  0, 
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(4.37) 


where  £N  =  e  +  d2(l  +  I)^1^2  >  *■  This  gives  us 

l!Pi(r)  -  P2(r) II  <  il$1  ~~  $2H  exP {-UnZnt}. 

In  addition,  we  have 


ll«i  -  U2II  <  (1  +  OilPi  ~  P2II  +  ||$i  -  *a|| i/m- 

Thus  combining  with  (4.23)  we  get 


|F(«i)  ~  -^(“2)1  <  d2  [  (l  +  0||Pi  -  P2II  +  ||$i  —  ^2 1|  • 


This  gives 


l|r$,(p„)  -  r*,(po)||  < 

[  HQiv5(-r)||J;(gNH;QwV)|  (l  +  0IIPl(T)  “  Pj(r)ll  +  11^1  "  $*11  VT ■ 

J  -oo 

By  applying  the  estimate  (4.37)  we  get 

l|r*i(p0)  -  r$i(p0)||  < 

rf2||$l  -  $l||  [  ||QN5(-7-)||£(gwH;gNV')[l  +  ~T^Nf2( 1  +  0e  l‘rji!Jt\dr. 

J  —  oo  t 

From  (4.27),  (4.35)  and  (4.36),  we  finally  get 

||T<Mp0)  -  r$2(Po)ll  <  rld2{  2e-1/V)v1+/i2  +  Pn1/2  OII*i  -  ^2 II 

=  L  ||$!  -  $2||, 


where  /'  is  given  in  (4.30). 

In  order  for  T  to  be  a  strict  contraction  we  must  have  L  <  1.  Let  us  choose 

L  =  d2e_1(2e“1/2/xN1/i  +  M*7*  H  <  \- 
Since  l'  <  l  and  n^2  <  pjy215  a  sufficient  condition  for  L  <  1/2  is 

p]v+ i  >  c_1d<,  Wlt^  =  2 d2  (  2e-1/2  +  / )  and  /  <  2. 


(4.38) 
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□ 

Let  us  now  establish  the  existence  of  the  inertial  manifold  for  (4.12). 

Theorem  4.1  Let  $  £  Hbi  with  0  <  /  <  1  and  b  be  given  in  (4.28).  Let  N*  be  given  and 
for  R  =  rz  +  b,  along  with  0  <  7  <  1  —  there  exist  constants  dz ,  d 4  depending  only  on 
and  H  such  that  the  following  three  conditions  are  satisfied: 

(i)  N  >  N*] 

(ii)  /$ii  > 

(iii)  VnIi  ~ 

Then  the  transform  T  has  a  fixed  point  $  £  Hbi.  The  manifold  .M  defined  by  the  graph  of 
$  is  an  inertial  manifold  for  (4-12). 

Proof:  From  Lemma  4.3  and  Lemma  4.7,  we  see  that  the  transform  T  maps  Hb  /  into 
itself  and  is  a  strict  contraction.  Thus  there  exists  a  unique  fixed  point  £  Hb ;  such  that 
r$*  =  <I>*  by  the  contraction  mapping  theorem. 

We  note  here  that  by  construction  the  manifold  M  is  invariant  to  the  action  of  the 
nonlinear  evolution  operator  Wm(t,0]  •).  To  see  this  we  first  write  the  relationship  T$*  = 
as 

$*(Po)  =  -  [  QNS(-T)QNF(u(r-,**,p0))dT,  (4.39) 

with  u(r,p0)  =  p(7-;$‘,p0)  +  $*(p(r;  $*, p0)). 

Let  p{t)  =  p(f;$*,p0)  be  the  solution  of  (4.15)  that  is  defined  for  all  t  £  R.  Then  we 
need  to  show  that  q(t )  =  $*(p(f))  is  a  solution  of  (4.14).  Let  us  consider  $*(p(<))  and 
using  the  fact  u(r,  p(t ))  =  u(r,  Wm(t,  0;  p0))  =  u(r  +  t ;  p0)  to  get 

$*(p(0)  =  -  f  QNS(-T)QNF{u{T,p{t)))dr 

J  —  OO 

=  -/  QnS{-t)Qn  F(u(t  +  t,p0))dr 

J  -00 
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By  setting  t  =  rj  —  t  and  from  (4.24),  (4.25)  one  has 


$*(p(t))=-f  QNS(t  -  r])QNF(u(r},p0))dri,  for  all  t  E  R. 

J  ~  oo 

By  making  use  of  Leibnitz  rule  formally,  we  find  that 

d$*(t)  r* 

-j.—  +  f-A Qn{  -  /  QNS(t  -  r})QivF(u(Tj,p0))drj}  +  QNF(u(t,p0 ))  =  0. 

at  J  —oo 

This  gives 

+  eAQ„**(p{t))  +  Q„F(p(t)  +  *>(*)))  =  0 

Clearly,  we  see  that  (p(t),  $*(p(t)))  is  a  solution  of  (4.13),  (4.14).  Hence  u(t)  is  a  solution 
of  (4.12)  with  q(t)  =  $*(p(t)).  This  proves  the  invariance  property  Wm(<,0;  >1)  C  M. 

We  shall  now  establish  exponential  attractive  property  of  the  manifold  M.  Let  us  recall 
a  theorem  on  the  squeezing  property  of  solution  orbits  [15]. 

Theorem  4.2  Let  Wm(t,0;u0 )  and  Wm(t,  0;»o)  be  two  solution  of  (4-12)  with  u0,  w0  E 
Br  C  V.  Then  there  exist  constants  depending  only  on  R,T,f, 

€,1/,*)  and  Cl  such  that  for  every  7  >  0  and  every  t  E  [0,  T],  we  have  either 

||C«M(WTO(t,0;uo)  -  Wm(t,0;v0)]||  <  7  0;  «0)  -  Wm(t,  0;  »0)]|| 


or 

||Wm(t,0;«0)  -  WmM;w0)U  <  exp(-C2  e  pN+it) ||u0  -  v0||  (4.40) 

for  every  N  >  1. 

We  will  now  show  that  this  property  implies  that  the  manifold  M  be  globally  exponen¬ 
tial  attracting.  For  convenience  we  choose  N *  such  that 

2C2  e  to 

This  gives 

||gw[Wm(<,0;u0)  -  Wm(t,0;«0)]||  <7||^[^m(<,0;u0)- Wm(<,0;Vo)l||  (4.41) 
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and 


||Wm(t,0;tio)  -  Wm(t,  0;  tr0)||  <  ^  ||u0  -  v0||  r4.42) 

for  t0  <  t  <  2t0. 

Let  us  denote  the  distance  between  any  point  w  in  the  absorbing  ball  BT3  in  V  and  the 
manifold  M  by 

dist(t»,  M)  ~  inf{jjtv  —  v||  :  v  E  At}. 

Let  v0  E  .M  so  that  v0  =  PNv 0  +  $(Pnv0)  and  that  ||u0  —  v0||  =  dist(«0,  .M).  Moreover,  it 
can  be  shown  that  ||PjvV0||  <  r3  for  v0  E  .M.  Hence  for  $  E  H^t,  it  follows  at  once  that 

H»o||  <  |!-Pn»o||  +  ||S(Pv*o)||  <  r3  +  b. 

hence  we  can  choose  R  =  r3  +  b  such  that  u0,v 0  E  Br.  Let  us  now  apply  the  squeezing 
property  stated  above.  We  shall  first  establish  the  attracting  property  in  to  <  r  <  2to 
using  either  one  of  conditions  (4. 41), (4. 42). 

Using  (4.42),  we  get 

dist(Wm(f,0;u0),  .M)  <  ||Wm(«,  0;  u0)  -  Wm{t,  0;  v0)|| 

<  i||«0  -  Vo||  =  X-  dist(«o,  X). 

Note  that  v0  will  always  stay  on  the  manifold  for  t0  <  t  <  2 10. 

Let  us  consider  a  point  «(<)  which  has  evolved  from  «0  E  Br.  That  is  u(t)  = 
Wm(t,0;«o)-  Then  there  is  a  point  on  the  manifold  X  which  can  be  represented  as 
v*  =  P^Wm(t,0;«0)  +  $(Pwiym(f,0;uo)).  It  follows  from  (4.41)  that 

dist(Wm(t,  0;  u0),  M)  <  \\Wm{t,  0;  u0)  -  v*|| 

<  ||<?NlUm(t,0;«0)  -  $(PNWm(t,0;u0))|| 

<  ||QwIUm(t,0;uo)  -  QNWm(t,0\  v0)|| 

+  ||$(PNWm(t,  0;  v0))  -  $(P;vWm(t,  0;  u0))  ||. 
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Since  M  is  invariant  to  the  action  of  Wm(t,  0;  •),  we  have  QNWm[t,  0;  Co)  '—$[PNWm(t,  0;  c0)). 
Note  that  we  also  have 


\MPNWm{t,  0;  Co))  -  *{PNWm{t,  0;  «0))  || 

<  l  \\PNWm{t,  0;  c0)  -  PNWm(t,0;  u0)||. 

From  this  it  follows  that 

dist(Wm(f,  0;  u0),  At)  <  (7  +  /)  \\PNWm{t,  0;  c0)  -  PNWm{t,  0;  u0)|| 

<  (7  +  l)\\Wm(t,  0;  c0)  —  Wm(t,  0;  u0)|| 

<  i||«0  -  »o||  =  ^  dist(u0,  M). 

Here  note  that  hn  <  M(v+i  and  7  +  /  <  1. 

From  the  semigroup  property  of  Wm(t,  0;  •),  we  can  deduce  that 

dist(VFm(f,0;u0),.M)  <  Q)  dist(u0,M). 

Suppose  we  write  t  —  nr  with  t0  <  r  <  2 10,  then 

dist(VFm(£,0;u0)>-M)  <  e~  5*oln2dist(u0,  >1). 

This  means  that 

dist(Wm(t,  0;  Uo),  M)  — 0  exponentially,  as  t  — *  00. 

4.4  Spectral  Growth  Rate 

Let  us  now  discuss  the  spectral  gap  condition  for  operator/!.  As  established 
previous  sections,  a  sufficient  condition  for  the  existence  of  inertial  manifold  is 

1/2  1/2  ^  -1  , 

for  each  e  >  0. 


□ 


in  the 


(443) 
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Let  us  consider  the  following  two  eigenvalue  problems: 

(— A)2^  =  i/tb 

. 

.  V’lan  =  0.  =  0 


(4.44) 


(-a)V  +  Vp  =  n<i> 

V  •  <t>  =  0 


(4.45) 


(  4>\an  =  0,  A0|an  =  0 

The  eigenvalues  Us  of  (4.44)  can  be  obtained  by  the  minimization  problem: 

.  (Au,Au)  .  f  |A«|2  y 

VN  =  mf  — - r—  =  inf  ueV, 

ueM^_,  (u,u)  ueMjJ.,  |u|2 

where  V  =  {u  6  H2( fl),u|an  =  0}  and  Ms- 1  =  span{Vq,  •  •  • ,  Similarly,  for 

eigenvalues  ns  of  (4.45),  we  minimize  the  Rayleigh  quotient  with  u  restricted  to  the  space 
V .  Since  V  is  a  subspace  of  V ,  by  First  Monotonicity  Principle  [22]  we  can  conclude 

Hn  >  vs- 

It  can  be  shown  that  the  growth  rate  of  the  eigenvalues  uN  of  (—A)2  for  an  arbitrary 
smooth  bounded  domain  f2  C  Rn  are  [4,17] 

..  .  ..  „  .  „ 


~ lfcr<  (bPv) 


as  N  — ►  oo. 


Here  V  is  the  volume  of  f2  and  Bn  is  the  volume  of  the  unit  ball  in  Rn.  Since  nN  >  vn, 


we  also  have 


c  Ni/n  as  N  -»  oo. 


Here  ns  is  eigenvalues  of  dissipation  operator  A  =  Ps(—A)2  and  n  is  the  space  dimension. 
In  particular,  for  space  dimension  n  =  2,  we  have  iiN  ~  c  iV2.  It  is  easy  to  see  that  for 


sufficient  large  N  we  get 


-  4’  ~  c‘/!- 


This  indicates  that  there  exists  an  inertial  manifold  only  for  e  >  —rrr  instead  of  each  e  >  0. 

c1/2 

Notice  that  if  eigenvalues  ns  ~  cN2+s  as  iV  — ►  oo,  then  the  large  spectral  gap  condition 
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in  (4.43)  is  satified  for  every  e  >  0.  This  will  be  true  if  A  =  PH(  —  A)2+S ,  8  >  0.  Hence  the 
dissipation  operator  A  associated  with  regularized  Navier-Stokes  equations  in  arbitrary 
two  dimensional  domains,  the  spectral  gap  condition  is  marginal.  However  the  spectral 
gap  condition  for  A  is  satisfied  for  the  following  case. 

Let  us  consider  problem  (4.45)  with  periodic  boundary  conditions.  Using  Fourier  series 
we  get  the  eigenvalues  as 

with  k  =  (ku  •••,*:„). 


In  particular,  for  n  —  2  we  have 


167T4 

~L*~ 


{k\  +  kl)\ 


Let  us  now  set  S( kl,k2)  =  k\  +  k\,  thus  M(*=i ,*2)  =  c‘S(kl  k2)  Wlth  c  =  167 r4/L4.  According  to 
magnitude  of  Siklki),  we  can  rewrite  S(ti  *2)  in  a  sequence  as  Si  <  S2  <  •  ■  ■  and  this  will 
establish  an  ordering  of  the  eigenvalues  as  Hi  <  <  •  ■  •,  with  =  c  Sir  Since  SN  is  sum 
of  two  squares,  we  can  apply  a  number  theoretic  result  by  Richards  [18]: 

1  -  M//2  =  c1/2(5JV+  i  -  SN)  >  \  log  SN  as  N  ->  oo. 

4 


Here  we  note  that  when  fl  is  periodic,  A  has  a  zero  eigenvalue.  Hence,  in  order  for  the 
earlier  theories  to  apply  we  should  modify  A  by  A  +  0,  (3  >  0. 


4.5  Regularity  of  Inertial  Varieties 

In  this  section  we  will  study  the  regularity  properties  of  inertial  manifolds  At.  We  will 
obtain  a  sufficient  condition  for  the  inertial  variety  to  be  C1.  We  will  prove  in  particular 
that  once  the  existence  of  inertial  variety  is  established  then  higher  dimensional  manifolds 
are  automatically  Cl. 

Definition  4.2  Let  Hh  be  the  subspace  of  and  satisfy: 

(0  ||D$  (p)  ||  £(Ph  V  ,Q  tfV)  —  b\i  bx  >  0,  VpePyvF, 

(ii)  HDfcfp,)  -  D${p2)\\c[pnV.<QnV]  <  Zillpn  -  p2]|,  Vpnp2  G  PNV. 
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Note  that  H is  complete  [12]  with  the  metric 


dist($,$')=  sup  |j$(p)  --  $'(p)||. 
pePNv 

In  the  Lemma  4.2  we  established  the  global  boundedness  and  Lipschitz  continuity  for 
the  modified  nonlinear  term  F(u).  We  will  now  show  similar  properties  for  the  Frechet 
derivative  [£>F](u). 

Lemma  4.8  For  $  £  if*,/  and  p1,p2  £  PnV ,  we  have 

(i)  \{DF}(u)h\  <  Ki||h||,  for  fixed  u£V\  V  h  e  V.  (4.46) 

(ii)  |[DF](«1)ft1-[DF](ti2)fc2| 

<  (1  +  l)(K2\\hi\\  c0|jh2||)  Hpj  —  p2||  +  Ki\\hi  —  /t2j|, 

VttLttj  €  V\  Vhuh2  £  V.  (4.47) 

Proof:  (i)  Note  that  the  Frechet  derivative  [Di2](u)  of  R{u)  is 

[DR\(u)  h  =  i/Aih  +  B(u,h)  +  B(h,u)  for  fixed  u  £  V.  (4.48) 

From  (4.18)  and  (4.19),  we  get 

I  [DiZ](u)  h\  <  u\Ayh\  +  |B(ti,fc)|  +  \B{h,u)\ 

<(!/  + 20811*11)11*11. 

This  gives 

|[2?FK«)*|<«ri(|H|)|[M](*)*| 

<  (i'  +  4c0ri)  ||h||,  VftG  F 

Thus  for  fixed  «  £  V,  ||  («) || £(v;/r)  £  K\  with  K\  =  v  +  4c0rx. 

(ii)  For  any  hi,h2  £  V,  We  have  from  (4.48), 

[ZLR](tii)  hi  -  [Z>i?](u2)  h2  =  vAi{hi  -  h2)  +  B(ui  -  u2,hi) 

+B(u2,  hi  —  h2 )  -f-  B[hi  —  h2  tij)  B[h2,ui  —  u2). 
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It  follows  from  (4.18),  (4.19)  that 


|  [DR^uJhi  -  [DR){u2)h2 1 

<  (1/  +  Co||u2||  +  c0||«i||)  ||fti  -  A2||  +  c0(|| Ax ||  +  (j/l2||)  |j«l  -  «2||. 

Now,  let  us  consider 

[DFKttO  h j  -  [DF}(u2)  h2  =  0ri(IM)[^K«i)  h,  -  h2 , 

for  the  following  three  cases.  First  suppose  that  ||«x||,  ||u2|(  >  2 rx,  then  we  immediately 
have  [tfFKux)*!  -  [DF}(u2)  k2  =  0.  If  ||mi}|  <  2rx  <  jjti2!|,  then  0ri(M)  =~  This  g'ives 

|  [jDF](u,)  hx  -  (L>F](u2)  h2\  <  i “ndlttill)  -  MIM)I  I  [M](tn)fcx| 

<  2r~1Ki\\ul  -  u2\\  ||fti||. 

For  the  case  ||«x||,  |ju2||  <  2 ru  we  have 

|[ DFjluJht-lDFlM  h2\ 

<|«.1(l|u,||)-*.,(||t.2||)|llDBl(«i)/.1| 

+#r,(||aI||)||Dfil(«,)(.1-(£>S](«2)'i!| 

<  [(2 r;‘K,  +  c„)||/.,||  +  c0||/i2||]  ||«i  -  «2||  +  K,\\h,  -  fc2||. 

Combining  these  results  gives, 

|  [DF](ux)  hi  —  [Z?F](u2)  h2| 

<  (i^2||/ii||  +  CoIlM)  llui  —  Uzll  FC 1 1 1 ^ i  —  ^2 1| )  (4.49) 

with  K2  =  2rf  lK\  +  c<> 

Writing  ux  —  u2  =  pi-p2  +  [$(px)-$(p2)]  for  $  €  Hbih  we  get  ||«i-u2|[  <  (l+OIIPi ~ P2II- 
Ilence,  for  any  hi,h2  £  V,  (4.49)  becomes 

|[Z?F](«i)A1-[PF](«2)/l2| 

<  (1  +  /)(/f2||Ai||  +  c0\\h2\\)  IIPi  -  P2II  +  Kill^i  -  M- 
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□ 

In  order  to  prove  the  existence  of  finite  dimensional  C 1  variety  At  using  contraction 
mapping  theorem,  we  establish  the 

Lemma  4.9  For  $  G  Hi /,  if  p]y2  j  —  pj]2  >  t~lK7,  pn+i  —  P N  >  e~2Ke  and  l  <  2  then  T 
maps  i/j1;  into  itself: 

T  :  Hi  —  fTj„ 

where  K7  =  max{d3,  lf3,  K*,}  and  Ke  =  KiK4(l  +  bi). 

Proof:  Let  us  first  show  that  D(T$)(p0)  is  globally  bounded  in  Z(PnV  ',Q  nV)  with 
respect  to  p0  G  PyvV\ 

Lemma  4.10  If  $  G  Hh  and  p0  G  PyvP,  then  [D(T$)](p0)  G  £[PnV',QnV).  Moreover, 

1  /2 

:'/  (5in  =  e(pAr+i  —  p-jv)  —  ffi(l  +  f>i)p^  >  0  then  we  have 

ll(O(r*)|(p0)||£(P„^,„K|  <  Mi, 

where  Mi  =  Ki(  1  +  6i)p^1/2  [e_1  +  (e  -  ^tvCi^)-1]  e_1/2  exP  '^2etL>  ant * 

£in  —  e  +  ffi(l  +  6i)pN1//2. 

Proof:  From  the  definition  of  T$(p0)  in  (4.26),  we  have  V/»  G  P/vP 

[•D(T$)](p0)/i  =  —  f°  QnS(-t)Qn[DF}(u)  o  (I  +  [D^}(p))  o{Dp}(p0)hdr. 

J  —  OO 

Note  that  from  Lemma  4.5,  we  have  (T$)(p0)G  This  means  (T<I>)(-)  is  a  nonlinear 

mapping  from PN V  into  Q^V.  Hence  the  Frechet  derivative,  [Z)(r$)](p0)  G  Z{PnV \QnV) 
for  fixed  p0  G  P;vV\  By  the  definition  of  [P(P$)](p0)  given  above,  such  result  can  be 
shown  as  follows:  Note  that  p(t)  =  p(t;#,p0)  implies  p(-)  :  P^V  t— >  PNV .  Thus  Frechet 
derivative  of  p(t)  at  p0  is  [Dp](p0)  G  £(P^V; PnV).  Now  note  that  $(•)  :  P^V  •-*  Q^V 
implies  [£>$](p)  G  Z(PnV\ QnV).  We  get 

[Dp\(Po)  +  [D$}(p)  o  [ Dp}(p0 )  G  Z{PNV-V). 
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Since  F(-)  maps  V  into  H,  we  obtain  [D /'](u)  1=1  C(V\H).  Consequently,  we  have 
Qx[DF](u)  £  £(V\QnH).  Combining  all  the  above  results  : 


QiV[£>i*’](«)  o  (I  +  [£>$](p))  o  \Dp}(p0)  £  £{PnV;QnH). 

Recalling  the  properties  of  the  holomorphic  semigroup  S(— r)  described  in  section  3.1,  we 
have  Q*S(-r)(-)  :  QNH  ~  QNV.  Hence  L>[T$](p0)  G  £(PNV;QNV). 

Let  us  now  set  ip  -  (I  +  [Z?$](p))  o  [Dp}(p0)h  £  V,  then  we  get 

ll[£(r*)](Pb)fc||  <  f  \\QNS(~T)\\ciQNH,QNV]\QN{DF}(n)rP\dT.  (4.50) 

J  —  OO 

First  we  need  to  estimate  |Qjv[£>.F](u)t/>|.  Since  ip  £  V,  we  have  from  (4.46)  and  the 
estimate  ||[£>$](p)||£(P„y;QNv)  <  bu 

\  [DF}(u)iP\  <  Kx !!(/  +  [U*](p)J  o  [Dp\(p0)h\\ 

<^i(l  +  MII  \Dp}(P0)h\\. 


Let  us  consider  g  =  \Dp\{pQ)h  £  P^V ,  for  h  £  P^V .  Taking  the  Frechet  derivative 
with  respect  to  p0  in  the  following  evolution  equation 

f  pt  +  eAp  +  PnF(u)  =  0 


{  P{t)  = 

we  obtain 

gt  +  eAg  +  PN[DF}{u)ip  =  0 
.  0(0)  =  h. 

Taking  inner  product  with  Ag  for  (4.51)  to  get 


^llffll2  +  <  Kx (1  +  bx)Hgll  \Agf. 


(4.51) 


Then,  by  using  the  fact  that  \Ag\2  <  Pat||0|[2  we  obtain 

\\q\\-^\\q\\  >  ~i tPN  -f  Kx{\  +  bi)fj,]j'}  ||g||2 

This  implies 


dt 


Iffll  +  ARvfijJItfll  >  0,  where  CiN  =  e  +  KX{1  +  bx)p 


-1/2 


N 
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Integrating  from  r  to  0  gives 


lff(r)ll  <  llff(°))l  exP {-»n£int} 

<  ll^ll  exp{-/ijvfi  t},  V  r  <  0. 


(4.52) 


Consequently,  we  have 


[DF}(u)xJ)\  <  K^l  +  6i)||ft||  exp{~nN^  t},  V  r  <  0. 


(4.53) 


By  applying  result  (4.53)  in  (4.50)  we  obtain 


||[B(r*)](Po)ft||  <  jt,(i  +  »,)l|ft|l /_°m  lie*s(->-)llt(«  NH,QnV)€-  »N*1NTdT 

<k,(i  +  MIWI  {(*!£./ 

—oo 

/  1  \ 1/2  r° 

+  ( -  /  |r|-l/2e-^flA,r(fr}> 

\2eeJ  J-i/2€iiN+i 


T\~1/2e~M,,(i"rdr}. 


Note  that  here  we  have  similar  estimates  as  in  (4.35),  (4.36).  i.e. 

Prill  [  1  +  eSlNTdr  <  e_1/?2  /*^Vi  (e  “  J_1  e" 


with  6lAf  =  €{hn+i  -  Hn)  ~  K\{^  +  t>i)VN,  oN  =  Hn/v-n+i  and 


1  \J/2  /■« 


n 


l/2f^/V  +  l 


|r|  x/2e  ,iN^lNTdT  <  e  1e  1^2  fJ-N+l  exP  2t  *  • 


We  then  have 


||[z?(r$)](p0)/»||  <M1||fc||,  vhePNv, 


with  Mi  =  Ki(l  +  bi)nN112  [e  1  +  (e  -  oNZ\N)  x]  e  1/2  exp  — 2^-.  Thus 


|j[I?(r$)](p0)|j,C(pJVV;QivV)  <  Mi. 


Let  us  now  prove  the  global  Lipschitz  continuity  of  map  [D(T$)}(-).  Note  that  here 
(D(T^)j(-)  defines  an  operator  valued  map  PnV  »  C(PnV',QnV)- 
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Lemma  4.11  If  6[N  =  tPN+i  ~  PnU in  +  6v)  >  0  then  for  $  G  and  P01,P02  €  PnV 
we  have 

II  |B(r*)|(Pol)  -  |-D(r$)|(p02)||i(p„V!Q„K)  <  /;  IIpoi  -  Poill. 

with 

l\  =  c3  1  +  K\[\  +  &i)e  1  nN  ^2] 

x  j^ArVi2  [e_1  +  (e  -  +  6v))_1]  e_1/2  exp  -  — |  • 

Proof:  Let  $  be  a  fixed  element  in  H^.  Let  pt(t)  and  p2(t )  be  two  different  solutions 
of  (4.15)  with  Pi(0)  =  p01  and  p2(0)  =  p02  respectively.  Let  us  first  consider  the  term 
[  [DF](ul)ip1  —  [DF](u2)il)2\,  from  (4.47)  we  have 

<  (i  +  OUW.II  +  «oll^j||)  Up.  -  p2|!  +  tf.ll*.  -  M,  (4.64) 

with  tt.  =  U  +  |-D$1(p,))  O  |Op.](Poi)ft  e  V  and  xp,  =  (I  +  (M](p,))  o  |Dpi](p0j)/>  €  V, 

V  h  €  PNV .  By  setting  g2  =  [OPiKPoi)^  e  P^V  and  g2  =  [Dp2\{pm)h  £  PNV ,  we 

immediately  have 

lhM<(i  +  *0IM  (4-55) 

and  ||02||  <  (l  +  Mll^ll-  (4-56) 

Since  -  02  =  (gx  ~  g2)  +  ([£>$](?!){/!  -  \D$]{p2)g2), 

1101  -  02ll  <  llffi  -  toll  +  II  P&|(Pi)?i  ~  [^Kft^all. 

where 

II  (M|(Pi)».  -  [r>«)(ft)»:ll 

<  ||  [D*](p,)(»,  -  »,)||  +  ||  ([B*|(P.)  -  [B#|(Pi))Sill 

<  tills.  -  Sill  +  (.l|Pi  -  Pill  llsill- 
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This  gives 


1101  -  02 II  <  (1  +  Mllfe  -  fell  +  *i||Pi  -  P2II  II  fe  II  ■  (4.57) 

By  applying  (4.55),  (4.56)  and  substituting  (4.57)  into  (4.54)  we  get 

|[DF](u1)01-[DF](ti2)02| 

<  Mfell  +  c2||fe||)  ||Pi  -  P2I!  +  KX{1  +  fci) ll^i  -  fell. 

where  cx  =  K2{  1  +  /)(1  +  6j),  c2  =  c0(l  +  /)(1  +  6j)  +  Kill.  From  (4.52),  note  that  we  also 
have  the  following  estimates: 

IlfeOOII  <  exp{-/xjv^lArr}  \\h\\,  V  r  <  0; 

IlfeOOII  <  exp{-pNtiNr}  ||*||.  V  r  <  0; 

IIPiOO  -  Pi(r)||  <  exp{-nN£NT}\\p0l  -  p02||  V  r  <  0. 


Now,  we  need  to  estimate  ||fe  — fe||.  Let  us  consider  the  following  two  evolutionary  problems 
for  ^(t)  and  g2(t)  respectively: 


9it  +  cA0i  +  Pn[DF](u1)iJ}1  =  0 

< 

.  tfi(o)  =  fc,  vfc  e  pKv, 


(4.58) 


and 


fe(  +  e44ff2  +  Ph[DF](u2)iI)2  =  0 
.  fe(o)  =  *»  Vh  €  PNV , 


(4.59) 


with  «i  =  Pi  +  $  (Pi)  and  «2  =  p2  4-  $(p2).  By  subtracting  (4.59)  from  (4.58),  we  get 


gt  +  eAg  +  Pn[DF}(ui)iJ>i  -  PN[DF](u2)il>2  =  0 

.  p(0  =  fe(0  -  fe(0>  ff(°)  =  °- 


(4.60) 


Taking  inner  product  of  Ag  with  (4.60)  to  get 


Ma  +  «M2 

<  |[L>F](u1)01-[L>F](u2)02|  \Ag\ 

<  hllPoi  -  P02II  ll*l|e_MwU,w+ew)‘  +  Kx[  1  +  61)||0||]  p]i2\\g\\, 
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which  implies 


^M\  +  /HvCijJull  >  -CzH^e  <*J»(«‘j»+€i*)‘|| pQl  -  p02 1|  \\h\\, 

with  C3  =  Ci  C2 .  This  leads  to 

llffi(0  -  fir2(r)ll  <  7^1/2HPoi  -  P02I!  exp{-fiN{CiN  +  &v)r}||/i||  V  r  <  0. 

Hence, 

l[£’f](«,)V>i-[£>F](»!)«I| 

<  CJ  [1  +  Jf,(l  +  6,)£“Vk1/2]  e“,‘''(£‘»+£'‘)l||p0,  -  P02 II  11*11- 

From  (4.61),  we  obtain 

II  [£>(r$)](p01)k  -  (D(r$)](p02)/>|| 

<  /  |lQp5(-T)||c(<3NH;QNv)ll^-f’](«l)^l  -'[r>P](U2)  ^,| 

J  —OO 

<  Cj  [l  +  Kt(l  +  6,)£"V«/I]  IIPoi  -  P02II  11*11 

J  -00 

where 

f  ll9*S(-r)||cw„Hi0„10e-'‘"<£'»+£''>'* 

y-00 

<  ^n+i  [c_1  +  (£  -  °n{Zin  +  6v))_1]  e-1/2  exp  gJV^1g£+  — 
with  =  e^N+i  -  Hn{Zin  +  Cn)  >  0-  Finally,  we  get 

|| [T>(r*)](p01)*  -  [z?(r$)](p02)/»||  <  /;i|p01  -  P02II IIAII,  v/*  g  p„v, 

with 

l'i  =  c3  [l  +  Ki(l  +  bi)t  Vjv^2] 

X  | P~nI\  [f-1  +  (f  -  °n{Zin  +  6v))_1]  e~1/2expg/^^1^+  —  |  . 

Thus, 

II  |B(T*)](p01)  -  |D(r*)|(p„,)||£(p„Vi«„v)  <  /;  llPo,  -  PoJ- 


□ 


Now  in  order  for  the  transform  T  to  map  Hlt  into  itself,  along  with  the  hypotheses  in 
Lemma  4.3  we  must  also  have  Mi  <  &i  and  l[  <  1%.  An  elementary  calculation  as  before 
show  that  a  sufficient  condition  for  Mi  <  b\  is 


1/2 

W+i 


1/2  .  - 
Vn  >  e 


lK, 


with  Kz  =  2K\(\  +  b{)bl  1 , 


(4.62) 


and  for  l[  <  li  is 


(  v]ll  i  -  h//2  >  c  1  ,  with  X5  =  max{X4, /fi(l  +  &i)  +  d2(l  +  /)} 

[  Vn+i  -  Vn  >  e~2K6,  with  K6  =  KiK4(l  +  bx),  K4  =  4c3/^1. 

Combining  above  results  (4.62),  (4.63)  with  the  hypotheses  in  Lemma  4.3  we  get 

Vn+i  -  Vn2  >  e~lK7,  with  K7  =  max{d3,  K3,  K5} 

< 

.  Vn+i  ~  Vn  >  t~2K$. 


(4.63) 


□ 


In  previous  section,  we  have  proven  that  T  is  a  strict  contraction  map  on  Htj.  Since 
Hh  is  a  closed  subspace  of  Hij,  this  immediately  implies  that  T  is  also  a  strict  contraction 
on  Hh-  Together  with  the  Lemmas  proved  above,  we  can  conclude  that  there  exists  a 
unique  fixed  point  such  that  by  the  contraction  mapping  theorem. 

Hence  we  have  the  following  theorem: 

Theorem  4.3  Let  4>  £  Hh  and  with  the  hypotheses  be  given  in  Theorem  4-1.  There 
exist  constants  d4,K6,K7  depending  only  on  l,li,buri,u,f  and  0  such  that  the  following 
conditions  are  satisfied: 


(i)  N  >  N*', 

(fi)  Villi  > 

(hi)  Villi  -  vT  >  t~lKr, 
(iv)  Vn+i  ~  Vn>  t~2K6. 
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Then  the  tranform  T  has  a  fixed  point  $  £  Hi i-  The  finite  dimensional  C 1  manifold  M 
defined  by  the  graph  of  $  is  an  inertial  manifold  for  (4-12). 


□ 


Let  us  now  compare  the  spectral  gap  conditions  for  Theorem  4.1  and  Theorem  4.3. 
Notice  that  for  a  fixed  e  >  0  there  exist  infinitely  many  choices  of  N  such  that  the  conditions 
in  Theorem  4.1  are  satisfied.  Let  us  assume  now  that  N0  is  the  smallest  number  for  which 
the  conditions  in  Theorem  4.1  are  satisfied.  Now,  for  the  existence  of  a  C1  inertial  manifold, 
additional  conditions  (iii)  and  (iv)  in  Theorem  4.3  are  needed.  Note  that  condition  (iv)  in 
Theorem  4.3  can  be  written  as 


Since  Pn+i  >  MNi  we  obtain 


1/2  1/2  s  1 

fiN+ 1  —  ~ 


K6 


e  V2c^JV+i. 


For  a  fixed  e  >  0  and  every  N  >  N0,  we  have 


i 

>  - 


Kq  i 


> 


1  K« 


« ' 2f  4b. 


Hence,  provided  that  N  >  N0,  condition  (iii)  and  (iv)  can  be  combined  into  one  condition. 


i.e. 

thili  ~  thl2  >  -K*  with  Kt  = 

€  Z(L\ 

This  means  for  N  >  No,  a  sufficient  condition  for  the  existence  of  a  C1  inertial  manifold  is 

..1/2  1/2  ^  -1  ts 

tlN+ 1  —  M/v  d.  e  /tg. 

Such  result  implies  that  the  higher  dimensional  inertial  manifolds  are  automatically  C1. 
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Appendix  A 

Time  Analyticity  of  Strong  Solution 


In  this  section  we  will  show  the  time  analyticity  of  the  strong  solution  using  a  method 
introduced  in  [21]  for  the  conventional  Navier-Stokes  equations.  Let  us  first  recall  an 
existence  theorem  provided  in  [15]  fur  the  evolutionary  equation  of  (4.11). 

Theorem  A.l  Let  n  <  6  and  u0  £  V,  /  £  L°°(0, T;H)  be  given.  Then  there  exists  a 
unique  solution  of  problem  (4-11),  which  satisfies  u  G  L2{0,T;  D(A))  Pi  C([0,  T\;V)  and 
u'  G  L2(0,  T;  H). 


□ 

Let  us  denote  C  as  the  complex  plane  and  Hq,Vq,D{Aq)  be  respectively  the  com- 
plexification  of  the  spaces  H,V,D(A). 

Theorem  A. 2  Let  f  be  an  Hc-holomorphic  function  on  D  ( a  neighborhood  of  the  positive 
real  axis  R+ )  and  is  bounded  from  D  into  He-  Then,  the  unique  solution  u  in  Theorem 
A.l  is  an  Hc-holomorphic  function  in  a  neighborhood  of  the  positive  real  axis. 

Proof  :  By  considering  projection  of  (4.11)  in  PnHq  we  get  the  complex  differential 
system 

— +  e.Attjv(f)  +  uPNAiUN{q)  +  PNB(uN($),uN($))  =  PNf  (A.l) 
uN(0)  =  PNu0.  (A.2) 
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The  above  system  is  equivalent  to  a  finite  system  of  complex  ordinary  differential  equations. 
Such  a  system  has  a  unique  solution  Ujv(c)  which  is  holomorphic  in  a  neighborhood  of  the 
origin  [3,  Theorem  8.1,  Chapter  l].  Let  us  now  obtain  an  estimate  on  the  size  of  this 
neighborhood  in  the  f-plane. 

By  taking  inner  product  of  (A.l)  with  Aun($),  we  get 

--^-|juAr(.se“?)jj2  +  e  cos  9\Au^(se'9)  j2  f  u  cos0jVA2«v('Se,#)|2 
2  as 

=  -Re{e,s (B(«Ar(se‘ff),  A«N(se'5))}  +  Re{e’#(/,  Auw(se,tf))}.  (A. 3) 


The  right  hand  side  of  above  equation  has  following  estimates: 


2  |(/,  Ati7v(5e,#))|  <  2  j/j  \AuN(se'e 
ecos  0  ,  , 

<  - \Aun |2  - 


■l/l2- 


2  1  ” 1  ecos  6 

The  estimate  for  trilinear  term  of  real  case  in  (2.10)  can  be  extended  to  the  complex  case 

|6(u, t7, w) |  <  c'1|«|1/2||«||1//2||v||1/,2|Aw|1/2|u7| 

Vug  Vq,  v  G  D(Ac)  w  G  He- 


This  gives  us 


2  |(B(uN(se,s),uiv(se^)),  AuN(se,s))| 

<  2c'1|uN(sei#)|1/2||uN(se^)||  |  AuN(seie)  |3/2 

<  ^I^Mse.*)!2  + 

2  (COSt?)° 

2  (cosc/J 

Thus  equation  (A. 3)  becomes 


—  ||u^(se*s)||2  +  ecos  0|Au^(se'^)|2  +  2i/ cos0|VA2u^(se",)|2 
ds 


e  cos  9 


(cos0): 


59 


sup  I!«iv(f)||2  <  ku  fci  =  fci(||u0||).  (A. 5) 

f6AD{0) 

Let  us  now  consider  the  passage  to  the  limit  of  the  solution  u^r(f).  From  the  estimate 
(A. 5)  and  the  compactness  of  the  embedding  Vic  C  He,  we  conclude  that  for  f  6  Ar>(0), 
{uw(f)}j£=1  belongs  to  a  compact  set  in  He-  Hence,  by  the  vector  valued  version  of 
the  Vitali  Convergence  Theorem  [7,  Theorem  8.2.1],  there  is  a  subsequence  {«jv,(<r)}“i  °f 
{uAf(?)}/v=1  that  uniformly  converges  to  JTc'holomorphic  function  u‘(f)  on  every  compact 
subset  of  A £>(0) .  Notice  that  the  restriction  of  ujy(f)  to  the  positive  real  axis  will  concides 
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with  the  Galerkin  solution  Ujv(t).  Hence,  it  is  also  true  that  the  restriction  of  ti*(<r)  to 
the  positive  real  axis  concides  with  the  unique  solution  u(t)  described  in  Theorem  A.l. 
Therefore,  u*(f)  can  be  viewed  as  a  holomorphic  extension  of  u(t)  to  the  region  Ab(0).  This 
implies  u'(f)  is  unique  as  well  in  the  holomorphic  region  Ap(0).  Due  to  the  uniqueness  of 
u'(f)  on  Ad(0),  the  entire  sequence  {un{$)}n=i  should  converge  to  u*(f)  in  Hq  uniformly 
on  every  compact  subset  of  A£)(0).  We  will  thus  denote  u*(f)  =  tx(<r). 

Let  us  now  consider  the  inequality  in  (A. 4)  for  0  <  to  <  a  and  s  >  t0.  Integrating  the 
differential  inequality  from  to  to  s  we  obtain 

||u^(se‘^)||2  <  const.,  with  t0  <  s  <  t0  +  a(cos 0)3. 

This  implies  that  U;v(f)  is  holomorphic  in  the  region  Ao(t0)  of  the  complex  plane  C,  and 

Ao(to)  =  D  n  A(t0)i  A (t0)  =  {c  =  se'e,t0  <  s  <  t0  +  a(cos0)3,  \0\  < 

£ 

Iterating  this  argument,  we  can  conclude  that  u(f)  is  an  JLc-holomorphic  function  in 
the  region  A  of  the  complex  plane  C,  with 

A  =  |J  Ar>(to)- 

<o>0 

That  is,  u(f)  is  an  L/c-holomorphic  function  in  a  strip  containing  the  positive  real  axis. 

□ 


61 


Bibliography 


[1]  R.  A.  Adams.  Sobolev  Spaces.  Academic  Press,  New  York,  1975. 

[2]  Garrett  Birkhoff  and  Gian-Carlo  Rota.  Ordinary  Differential  Equations.  Blaisdell 
Publishing  Company,  Waltham,  Mass.,  1969. 

[3]  E.  A.  Coddington  and  N.  Levinson.  Theory  of  Ordinary  Differential  Equations. 
McGraw-Hill,  New  York,  1955. 

[4]  R.  Courant  and  D.  Hilbert.  Methods  of  Mathematical  Physics,  volume  1.  Interscience 
Publishers,  Inc.,  New  York,  1953. 

[5]  J.  Dieudonne.  Foundations  of  Modern  Analysis.  Academic  Press,  New  York,  1969. 

[6]  C.  Foias,  G.  R.  Sell,  and  R.  Temam.  Inertial  manifolds  for  nonlinear  evolutionary 
equations.  J.  of  Differential  Equations,  73:309-353,  1988. 

[7]  E.  Hille.  Method  in  Classical  and  Functional  Analysis.  Addison-Wesley  Publishing 
Company,  1972. 

[8]  E.  Hopf.  A  mathematical  example  displaying  the  features  of  turbulence.  Comm.  Pure 
Appl.  Math.,  1:303-322,  1948. 

[9]  T.  Kato  and  H.  Fujita.  On  the  nonstationary  Navier-Stokes  system.  Rend.  Sem. 
Univ.  Padova ,  32:243-260,  1962. 

[10]  O.  A.  Ladyzhenskaya.  A  dynamical  system  generated  by  the  Navier-Stokes  equations. 
J.  Soviet  Math.,  3(3):458-479,  1975. 


[  1 1  ]  O.  A.  Ladyzhenskaya.  Mathematical  analysis  of  Navier-Stokes  equations  for  incom¬ 
pressible  liquids.  Annual  Rev.  of  Fluid  Mech.,  7:249-272,  1975. 

[12]  0.  E.  Lanford.  Bifurcation  of  periodic  solutions  into  invariant  tori ,  volume  322  of 
Lecture  Notes  in  Mathematics ,  pages  159-192.  Springer-Verlag,  Berlin-Heidelberg, 
New  York,  1973. 

[13]  J.  L.  Lions.  Espaces  d’interpolation  et  domaines  de  pussances  fractionnaires 
d’operateurs.  J.  Math.  Soc.  Japan ,  14(2):233-241,  1962. 

[14]  J.  L.  Lions  and  E.  Magenes.  Nonhomogeneous  Boundary  -Value  Problems  and  Appli¬ 
cations.  Springer,  Berlin,  1972. 

[15]  Y.-R.  Ou.  Analysis  of  regularized  Navier-Stokes  equations.  PhD  thesis,  University  of 
Southern  California,  Los  Angeles,  October  1988. 

[16]  A.  Pazy.  Semigroups  of  Linear  Operators  and  Applications  to  Partial  Differential 
Equations.  Applied  Mathematical  Sciences,  44.  Springer-Verlag,  Berlin,  Heidelberg, 
New  York,  1983. 

[17]  A.  Pleijel.  On  the  eigenvalues  and  eigenfunctions  of  elastic  plates.  Comm.  Pure  Appl. 
Math.,  3:1-10,  1950. 

[18]  Ian  Richards.  On  the  gap  between  numbers  which  are  the  sum  of  two  squares.  Adv. 
Math.,  46:1-2,  1982. 

[19]  G.  R.  Sell.  Topological  Dynamics  and  Ordinary  Differential  Equations.  Van  Nostrand- 
Reinhold,  London,  1971. 

[20]  S.  S.  Sritharan.  The  invariant  manifold  theory  for  hydrodynamic  transition,  submit¬ 
ted  to  Acta  Math.,  June  1988. 

[21]  R.  Temam.  Navier-Stokes  Equations  and  nonlinear  functional  analysis.  CBMS-NSF 
Regional  Conference  Series  in  Applied  Mathematics.  SIAM,  Philadelphia,  1983. 


63 


[22]  Hans  F.  Weinberger.  Variational  Methods  for  Eigenvalue  Approximation.  SIAM, 
Philadelphia,  1974. 

[23]  S.  D.  Zaidman.  Abstract  Differential  Equations.  Pitman  Advanced  Publishing  Pro¬ 
gram,  San  Francisco,  London,  Melbourne,  1979. 


64 


NASA 


Report  Documentation  Page 


1.  Report  No  2  Government  Accession  No 

NASA  CR- 18 1801 

ICASE  Report  No.  89-14 

4  Title  and  Subtitle 


I  3  Recipient's  Catalog  No 


i 

—  -4 - - 

;  5.  Report  Date 


ANALYSIS  OF  REGULARIZED  NAVIER-STOKES 
EQUATIONS  -  II 


I  February  1989 

1  6  Performing  Organization  Code 


i 


I 

1 


i 

i 


Authorise 


8  Performing  Organization  Report  No 


Yuh-Roung  Ou 
S.  S.  Sritharan 


!  89-14 

I  10  Work  Unit  No 


9  Performing  Organization  Name  and  Address 


Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 
PstRAtpAgeftc v 


|  National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  VA  23665-5225 

j  15.  Supplementary  Notes 

i 

Langley  Technical  Monitor:  Pacific  J. 

,  Richard  W.  Barnweil 


1  505-90-21-01 

11  Contract  or  Grant  No 

j  NAS 1-1LoU5 

13.  Type  of  Report  and  Per-d  Covered 

Contractor  Report 
14  Sponsoring  Agency  Co.'e 


of  Mathematics 


I  Final  Report 
!  16  Abstract 

j 

I  A  practically  important  regularization  of  the  Navier-Stokes  equations  have 

been  analyzed.  As  a  continuation  of  the  previous  work,  wa — study  in  this  paper 
the  structure  of  the  attractors  characterizing  the  solutions.  Local  as  well  as 
j  global  invariant  manifolds  have  been  found.  Regularity  properties  of  these 
manifolds  are  analyzed. 


j 


j  17  Key  Words  (Suggested  by  Author(s)) 

18  Distribution  Statement 

1  ! 

[  •’“Navier-Stokes  equations, ' 
invariant  manifold’,  global 

attractor 

/ 

64  -  Numerical  Analysis 

34  -  Fluid  Mechanics  & 

Heat  Transfer 
Unclassified  -  Unlimited 

19  Security  Classif  (of  this  report! 

20  Security  Classif  (of  this  pagel 

21  No  of  pages 

22  Price 

Unclassi fied 

Unclassi fied 

68 

A0  4 

NASA  FORM  1626  OCT  86 


N  ASA-Langley,  1989 


